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Abstract 

We review in detail the construction of all stable static fermion bags in the 
1 + 1 dimensional Gross-Neveu model with N flavors of Dirac fermions, in the large N 
limit. In addition to the well known kink and topologically trivial solitons (which corre- 
spond, respectively, to the spinor and antisymmetric tensor representations of 0(2N)), 
there are also threshold bound states of a kink and a topologically trivial soliton: the 
heavier topological solitons (HTS). The mass of any of these newly discovered HTS's is 
the sum of masses of its solitonic constituents, and it corresponds to the tensor product 
of their 0(2N) representations. Thus, it is marginally stable (at least in the large N 
limit). Furthermore, its mass is independent of the distance between the centers of its 
constituents, which serves as a flat collective coordinate, or a modulus. There are no 
additional stable static solitons in the Gross-Neveu model. We provide detailed deriva- 
tion of the profiles, masses and fermion number contents of these static solitons. For 
pedagogical clarity, and in order for this paper to be self-contained, we also included de- 
tailed appendices on super symmetric quantum mechanics and on reflectionless potentials 
in one spatial dimension, which are intimately related with the theory of static fermion 
bags. In particular, we present a novel simple explicit formula for the diagonal resolvent 
of a reflectionless Schrodinger operator with an arbitrary number of bound states. In 
additional appendices we summarize the relevant group representation theoretic facts, 
and also provide a simple calculation of the mass of the kinks. 
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1 Introduction 



Many years ago, Dashen, Hasslacher and Neveu (DHN) [1], and following them Shei 
[2], used inverse scattering analysis [3] to find static fermion-bag [4, 5] soliton solutions 
to the large-iV saddle point equations of the Gross-Neveu (GN) [6] and of the 1 + 1 
dimensional, multi-flavor Nambu-Jona-Lasinio (NJL) [7] models. In the GN model, 
with its discrete chiral symmetry, a topological soliton, the so called Callan-Coleman- 
Gross-Zee (CCGZ) kink [8], was discovered prior to the work of DHN. In this paper 
we will concentrate exclusively on the GN model. 

One version of writing the action of the 1 + 1 dimensional GN model is 

S = J d 2 x\^^ a (i$-a)^ a -^ 2 a 2 Y (1.1) 

where the ip a (a = 1, . . . , N) are N flavors of massless Dirac fermions, with Yukawa 
coupling to the scalar 1 auxiliary field a(t,x). 

An obvious symmetry of (1.1) with its N Dirac spinors is U(N). Actually, (1.1) 
is symmetric under the larger group 0(2N). It is easy to see thisfl] in a concrete 
representation for 7 matrices, which we choose as the Majorana representation 

7° = &2 , 7 1 = i&3 and 7 5 = — 7°7 1 = U\ . (1.2) 

(Henceforth, in this paper we will use this representation for 7 matrices in all explicit 
calculations.) 

In order to expose the 0(2N) symmetry, we write each Dirac spinor as 

lpa = <fra + iXa (1.3) 

where <f) a and \a are hermitean two-component Majorana spinors. In terms of (1.2) 
and (1.3) we can write the lagrangian in (1.1) (up to surface terms) as 2 

N ■ a 2 



o=l 



1 The fermion bag solitons in these models arise, as is well known, at the level of the effective 
action, after integrating the fermions out, and not at the level of the action (1.1). 

2 Here we use the common notation where an overdot stands for a time derivative and a prime 
stands for a spatial derivative. 
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which is hermitean and non-vanishing, because all spinors are Grassmann valued. 
Evidently, (1.4) is invariant under orthogonal transformations of the 2N Majorana 
spinors <p a and Xa- 

The fact that the symmetry group of (1.1) is 0(2N) rather than U(N), indicates 
that (1.1) is invariant against charge-conjugation. In the representation (1.2), charge- 
conjugation is realized simply by complex conjugation of the spinor 3 



with energy u, then ip*(x) = e lbJt u*{x) is an energy eigenstate of (1.6), with energy 
— lu. Therefore, the GN model (1.1) is invariant against charge conjugation, and 
energy eigenstates of (1.6) come in ±u pairs. 

The remarkable discovery DHN made was that all self-consistent static bag config- 
urations in the GN model were reflectionless. That is, the static cr(x)'s that solve the 
saddle point equations of the GN model are such that the reflection coefficient of the 
Dirac equation (1.6) vanishes identically for all momenta. In other words, a fermion 
wave packet impinging on one side of the potential a(x) will be totally transmitted 
(up to phase shifts, of course). 

We note in passing that besides their role in soliton theory [3], reflectionless po- 
tentials appear in other diverse areas of theoretical physics [9, 10, 11, 12]. For a 
review, which discusses reflectionless potentials (among other things) in the context 
of supersymmetric quantum mechanics, see [13]. 

Since the works of DHN and of Shei, these fermion bags were discussed in the 
literature several other times, using alternative methods [14]. For a recent review on 
these and related matters (with an emphasis on the relativistic Hartree-Fock approx- 
imation), see [15]. 

3 The matrices 7°,7 1 are pure imaginary. Thus by taking the complex conjugate of the Dirac 
equation [i0 — eA — c(a;)] ^ = 0we conclude that [i0 + eA — &{x)\ ip* = 0. 



ij c (x) = i>*(x) . 



(1.5) 



Thus, if ip — e lwt u(x) is an eigenstate of the Dirac equation 



i0 — a(x) tp — , 



(1.6) 
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In many of these treatments, one solves the variational, saddle point equations 
by performing mode summations over energies and phase shifts. An alternative to 
such summations is to solve the saddle point equations by manipulating the resolvent 
of the Dirac operator as a whole, with the help of basic tools of Sturm-Liouville 
operator theory. The resolvent of the Dirac operator takes care of mode summation 
automatically. 

Some time ago, such an alternative to mode summation techniques was devel- 
oped [16, 17], which was based on the Gel'fand-Dikii (GD) identity [18] (an identity 
obeyed by the diagonal resolvent of one-dimensional Schrodinger operators, see e.g. 
(A. 33)) 4 , to study fermion bags in the GN model [16] as well as other problems [17]. 
In a nut-shell, the method of [16, 17] is based on the fact that in certain models, 
the explicit form of the saddle-point equations for one-dimensional (or quasi one- 
dimensional) static, space-dependent field condensates, suggests certain parameter 
dependent ansatzes for the diagonal resolvent of the Dirac (or Klein-Gordon) opera- 
tor. This explicit construction of the diagonal resolvent is based on the GD identity as 
well as on simple dimensional analysis. All subsequent manipulations with these ex- 
pressions involve the space dependent condensates directly, and given such an ansatz 
for the diagonal resolvent, one can construct a static space-dependent solution of the 
saddle-point equations in a straightforward manner, bypassing the need to work with 
the scattering data and the so called trace identities that relate them to the space 
dependent condensates. 

That method was applied in [16, 19, 20] to study static fermion bags in the GN 
and NJL models and reproduced the static bag results of DHN and of Shei. It was 
also used in [21] to make approximate variational calculations of static fermion bags 
in the massive GN model [21]. Similar ideas were later used in [22] to calculate the 
free energy of inhomogeneous superconductors. 

In this paper we study the entire spectrum of stable static fermion bags in the GN 
model in the large N limit, thus extending the "static part" of the work[l] of DHN. 
4 For a simple derivation of the GD identity, see [17, 19]. 
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We find that in addition to the well known kink and topologically trivial solitons, 
which correspond, respectively, to the spinor and antisymmetric tensor representa- 
tions of 0(2N), the GN model bears also threshold bound states of a kink and a 
topologically trivial soliton - the heavier topological solitons (HTS)[23]. The mass of 
any of these newly discovered HTS's is the sum of masses of its solitonic constituents, 
and it corresponds to the tensor product of their 0(2N) representations. Thus it is 
marginally stable (at least in the large N limit). Furthermore, the mass of an HTS is 
independent of the distance between the centers of its constituents, which thus serves 
as a flat collective coordinate, or a modulus. There are no additional stable static 
solitons in the Gross-Neveu model. 

The rest of the paper is organized as follows: In Section 2 we recall some basic 
facts about the GN model and its dynamics. In particular, we define the effective 
action S e ff[a], obtained from (1.1) after integrating the fermions out, and derive the 
generic saddle-point equation for a(x,t). 

In Section 3 we focus on the simpler problem of finding extremal static cr(x) 
configurations, and as we mentioned earlier, DHN have already shownfl] that the 
extremal static configurations are necessarily reflectionless. Our analysis employs 
extensively one dimensional supersymmetric quantum mechanics and the theory of 
reflectionless potentials. We apply these formalisms to construct a generic static 
reflectionless background a(x) (subjected to the appropriate boundary conditions at 
x = ±oo), which supports a given number of bound states of the Dirac equation 
(1.6) at some given energies, and the diagonal resolvent of the Dirac operator in 
(1.6) associated with it. We then use this explicit resolvent to calculate the energy 
functional associated with this reflectionless background, as well as its fermion number 
spectrum. A reflectionless background a(x) depends on a finite set of real parameters, 
and the effective action evaluated at such a a(x) configuration is an ordinary function 
of these parameters. By solving the ordinary extremum problem for this function we 
identify all static extremal fermion bags in the GN model. 

After listing all extremal static a(x) configurations in Section 3, we identify, in 
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Section 4, the stable configurations, namely, those extremal bags which are protected 
by the conservation of topological charge and 0(2N) quantum numbers against de- 
caying into lighter fermion bags. These are the well known CCGZ kink and DHN 
solitons, and in addition, our newly discovered marginally stable HTS. 

Some general background material, as well as many technical details are relegated 
to the Appendices. In Appendix A we discuss (in some detail) the general properties 
of the resolvent of the Dirac operator in a static a(x) background and its relation 
with supersymmetric quantum mechanics. As is well known, the Dirac equation in 
any such background is equivalent to a pair of two isospectral Schrodinger equations 
in one dimension, namely, a realization of one-dimensional supersymmetric quantum 
mechanics. We use this underlying supersymmetry to express all four entries of 
the space-diagonal Dirac resolvent (i.e., the resolvent evaluated at coincident spatial 
coordinates) in terms of a single function: the diagonal resolvent of one of the two 
isospectral Schrodinger operators. Furthermore, using the underlying supersymmetry 
of the Dirac equation, we prove that the expectation value of the spatial component 
of the fermion current operator in the static a(x) background vanishes identically, 
as should be expected on physical grounds. We also prove an identity relating the 
expectation values of the fermion density and the pseudoscalar density operators 
in the <j(x) background, and show that it has a simple interpretation in terms of 
bosonization. This appendix follows, in part, our recent discussion in [20]. 

In Appendix B we summarize some useful properties of reflect ionless Schrodinger 
hamiltonians and their resolvents. Following that summary of results from the lit- 
erature, we use them to derive an explicit simple representation for the diagonal 
resolvent of a reflectionless Schrodinger operator with a prescribed number of bound 
states, which we have not encountered in the literature. We then use the formalism 
reviewed in that appendix to derive an explicit formula for a cr(x) with a prescribed 
set of bound states. We also work out explicitly the cases of backgrounds with a 
single bound state and with two bound states, which are the cases most relevant to 
this work. 
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Appendix C contains some technical details concerning the derivation in Section 
2 of the saddle point equation in the static case. 

In Appendix D we discuss the quantization of a multiplet of 2N Majorana fields in 
a given static a(x) background and the related 0(2N) multiplet structure of bound 
states of fermions trapped in that background. We show there that each bound state 
in the corresponding Dirac-Majorana equation, at non-zero frequency, gives rise to 
an irreducible factor of an 0(2N) antisymmetric tensor, and that the zero mode, 
which exists if and only if a(x) has non-trivial topology, gives rise to a single factor 
of the spinor representation. Also included in that appendix are the necessary facts 
concerning the spinorial and antisymmetric tensor representations of 0(2N). 

In Appendix E, we give a simple derivation of the mass formula of the CCGZ 
kink. Finally, in Appendix F, we present an alternative proof of Eq. (4.12), which is 
the basis for the complete classification of stable static solitons. 
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2 Dynamics of the Gross-Neveu Model 



We now recall some basic facts about the dynamics of the GN model. The partition 
function associated with (1.1) is 5 

Z = J VaVipViP expi J d 2 x j$ (i$ - o) $ - ^ a 2 j (2.1) 

Integrating over the grassmannian variables leads to Z — J Do exp{iS e ff[a}} where 
the bare effective action is 

Seff[<j] = -^J d 2 xa 2 -iNTr\og(i0-a) (2.2) 

and the trace is taken over both functional and Dirac indices. 

The theory (2.2) has been studied in the limit N — > oo with Ng 2 held fixed[6]. In 
this limit (2.1) is governed by saddle points of (2.2) and the small fluctuations around 
them. (In this paper, as in [1], we will consider only the leading term in the 1/N 
expansion, and thus will not compute the effect of the fluctuations around the saddle 
points.) The most general saddle point condition reads 
SS eff a(x,t) 



+ iN tr 



= 0. (2.3) 



5a(x,t) g 2 

In particular, the non-perturbative vacuum of (1.1) is governed by the simplest 
large N saddle points of the path integral associated with it, where the composite 
scalar operator ipip develops a space-time independent expectation value. The relevant 
object to discuss in this context is the effective potential V e ff associated with (1.1), 
namely, the value of —S e ff for space-time independent a configurations per unit time 
per unit length. V e ff(o~) has two degenerate (absolute) minima at a — ±m ^ 0, 
where m is fixed by the (bare) gap equation [6] 

-m + ^a 2 tr / d k „ 1 — = (2.4) 
y J (2n) 2 U-m K J 



which yields the fermion dynamical mass 



m = Ae N a' 2 w . (2.5) 



3 From this point to the end of the paper flavor indices are usually suppressed. 
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Here A is an ultraviolet cutoff. The mass m must be a renormalization group invariant. 
Thus, the model is asymptotically free. We can get rid of the cutoff at the price of 
introducing an arbitrary renormalization scale \i. The renormalized coupling g R (/x) 
and the cutoff dependent bare coupling are then related through 



l- 



Ae = N9 r m , (2.6) 

in a convention where Ng R (m) = n. Trading the dimensionless coupling g R for 
the dynamical mass scale m represents the well known phenomenon of dimensional 
transmutation. 

In terms of m, we can write down the renormalized effective potential (in the large 
N limit) as 

V; //M = ^lo g ^. (2.7) 

It is evidently symmetric under a — > —a, which generates (together with ip — > 75^ in 
(1.1)) the discrete (or ^ 2 ) chiral symmetry of the GN model. 

This discrete symmetry is dynamically broken by the non-perturbative vacuum, 
and thus there is a kink solution [1, 8, 16], the CCGZ kink mentioned above, <r(x) = 
mtanh(mx), interpolating between the two degenerate minima a = ±m of (2.7) at 
x = ±00. Therefore, topology insures the stability of these kinks. 

The CCGZ kink is one example of non-trivial excitations of the vacuum, which are 
described semiclassically by large iV saddle points of the path integral over (1.1) at 
which a develops space dependent, or even space-time dependent expectation values. 
These expectation values are the space-time dependent solution of (2.3), and have 
analogs in other field theories [12, 24]. Saddle points of this type are important also 
in discussing the large order behavior [25, 26] of the expansion of the path integral 
over (1.1). 

These saddle points describe sectors of (1.1) that include scattering states of the 
(dynamically massive) fermions in (1.1), as well as a rich collection of bound states 
thereof. These bound states result from the strong infrared interactions, which polar- 
ize the vacuum inhomogeneously, causing the composite scalar ipif) field to form finite 
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action space-time dependent condensates. We may regard these condensates as one 
dimensional fermion bags [4, 5] that trap the original fermions ( "quarks" ) into stable 
finite action extended entities ("hadrons"). 

Finding explicit space-time dependent solutions of (2.3) is a very difficult problem. 
In [1], DHN managed to guess such a set of solutions which oscillate in time (and 
thus are not merely boosts of static bags). Finding a systematic method to generate 
space-time dependent solutions of equations like (2.3) is, of course, one of the basic 
goals of quantum field theory. 

2.1 Static Inhomogeneous a(x) Backgrounds 

In this paper we focus on the easier problem of solving (2.3) for static inhomogeneous 
condensates. 

Static solutions of (2.3) are subjected to certain spatial asymptotic boundary 
conditions. For the usual physical reasons, we set boundary conditions on our static 
background fields such that a(x) starts from one of its vacuum expectation values 
a = ±m at x = — oo, wanders along the a axis, and then relaxes back to one of its 
vacuum expectation values at x = +oo: 

a — ► ±m , a' — ► 0. (2.8) 

These four possible asymptotic behaviors determine the topological nature of the 
condensate. The two cases in which the asymptotic values are of opposite signs are 
topologically non trivial. The CCGZ kink (where, by definition, cr(oo) = m ) and 
antikink (er(oo) = —m) are two such examples. The other two cases, i.e., the cases 
with asymptotic values of equal signs, are topologically trivial. 

The topological charge associated with these boundary conditions is simply 

1 1 °° 

q = 2^ (a(00) ~ <J(_00)) = 2^ / d ' Mx) ' (2 ' 9) 

— oo 

and can take on values q — ±1 for the topologically non trivial configurations, and 
q = for the topologically trivial ones. 
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As typical of solitonic configurations, we expect that a(x) tends to its asymptotic 
boundary values (2.8) at an exponential rate which is determined, essentially, by the 
mass gap m of the model. It is in such static backgrounds cr(x) that we have to invert 
the Dirac operator and calculate its resolvent in seeking static solutions of (2.3). 

In addition to the CCGZ kink, there are also topologically trivial inhomogeneous 
condensates. These topologically trivial condensates are stable because of the binding 
energy released by the trapped fermions, and therefore cannot form without such 
binding. Thus, they are stable due to dynamics. In contrast, the stability of the 
CCGZ kink is guaranteed by topology already. It can form without binding fermions. 

Note in passing, that the 1 + 1 dimensional NJL model, with its continuous symme- 
try, does not have a topologically stable soliton solution. Thus, unlike the GN model, 
the solitons arising in the NJL model can be stabilized only by binding fermions 
[2, 19]. This description agrees with the general physical picture drawn in [27]. 

2.1.1 The Energy of a Static Configuration and Supersymmetry 

In order to study the extremum condition (2.3) on S e g around a static space depen- 
dent background cr(x), we need to invert the Dirac operator 

D = ivy + i^ l d x — a(x) . (2.10) 

(We have naturally transformed %$ — a(x) to the cu plane, since a{x) is static.) In 
particular, we have to find the diagonal resolvent of (2.10) in that background. 

There is an intimate connection between the spectral theory of the Dirac oper- 
ator (2.10) and one dimensional supersymmetric quantum mechanics [13, 16]. This 
connection is explained in detail Appendix A. In particular, the topological charge q 
(2.9) coincides with the Witten index of the one dimensional supersymmetric quan- 
tum mechanics (see Eq. (A. 22)). 

The desired diagonal resolvent (x |i-D _1 |:r) is defined in (A. 29) in Appendix A. It 
is shown there, that due to the underlying supersymmetric quantum mechanics, its 
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four entries 

A(x) B(x) \ 
C(x) D(x) J 

are completely determined, in a simple way, in terms of any one of the off-diagonal 



(x | — iD \x 



(2.11) 



entries, B or C . 

For example, in terms of B, one finds (see Eqs. (A. 30), (A. 35) and (A. 38)) 



A(x) = D(x) 



j \d x + 2a{x)]B{x) 



-u 2 C(x) = =-B" + aB' + (a' + uj 2 )B . 



(2.12) 



A quick way to introduce this underlying supersymmetry is the following: From 
the elementary identity 

75(^-^)75 = -(i# + <x) (2.13) 

it follows that Trlog(i($ — a) = |Trlog — (i0 — a)(i$ + a) . Thus, an alternative 
representation of (2.2) is [16] 



5 e//M = -7^ / ^^-^Trlog^ + a 2 -^^) . 



2g* J - ~ - '2 

If, in addition, a is time independent, (2.14) may be further simplified to 



(2.14) 



SeffW] 



T 



dxa — i 



NT rdco 



Tr\og(H b -u 2 ) + Tr\og(H c -u 2 ) , (2.15) 



2g 2 J "" ' 2 J 2tx 

where T = / dt is a large temporal infrared cutoff, and 

H b = -dl + a 2 -a' and H c = -dl + a 2 + a' 



(2.16) 



are the pair of isospectral Schrodinger operators defined in (A. 6). In other words, 
they have the same spectrum of bound state energies, save, possibly, the existence 
of a bound state at u 2 = 0, which, if it exists, appears in the spectrum of only one 
of the operators. As explained in subsection A. 1.1 of Appendix A, a bound state at 
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uj 2 = exists only in topologically non-trivial a(x) backgrounds. Due to the zero- 
mode mismatch of the spectra of H b and H C) we will keep both Tr \og(H b — uj 2 ) and 
Tr \og(H c — uj 2 ) explicitly in our formulas for a while. We will take advantage of the 
isospectrality of H b and H c in the positive part of the spectrum at the appropriate 
places. 

The operators H b and H c may be identified as the hamiltonians of the bosonic 
sector and of the fermionic sector of a one dimensional supersymmetric quantum 
mechanical system. 

According to (A. 31), the off-diagonal entries B(x) and C(x) in (2.11) are essen- 
tially the diagonal resolvents of H b and H c : 

B(x) | 1 , . 

C{X) ~ <*l^l*>- (2-17) 



uj H c — uj 



From (2.15) we may express the energy functional £[<7(x)] of a static configuration 
a(x) as 



S[a(x)] = -^^ = ^- 2 Jdxa 2 + i ^ [Trlog(^-c 2 )+Trlog(/J c -^ 

9 71 (2-18) 

This expression is divergent. We regulate it, as usual, by subtracting from it the 
divergent contribution of the vacuum configuration a 2 = m 2 and by imposing a UV 
cutoff A on uj. Thus, the regulated (bare) energy functional associated with a(x) is 



tr°[<r{x)\ = ^- 2 jdx(a 2 -m 2 ) +l ^ J ^-{Tr\og(H b -uj 2 )-Tr\og(H VAC -uj 2 ) 

M<A 

+ 4 / % [ Tr l0g( ^ c - ^ - Tr Xo ^ H vac - ^ 2 )] , (2-19) 



where 



H VAC = ~dl + m 2 (2.20) 
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is the hamiltonian corresponding to the vacuum configuration. We are not done 
yet, as the integrals over u in (2.19) still diverge logarithmically with the UV cutoff 
A. However, as will be clear from the explicit calculations in the next section, this 
divergence is canceled by the logarithmic A dependence of the bare coupling g 2 , as 
determined by (2.5) and (2.6). 

The renormalized quantity £[<7(x)] thus defined in (2.19) is the mass of the static 
fermion bag. 
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3 Extremal Static Fermion Bags 

Now that we have written the energy-functional (2.19) of a static configuration, or 
a fermion bag, our next step is to identify those fermion bags on which (2.19) is 
extremal. 

The energy functional (2.19) is, in principle, a complicated and generally unknown 
functional of a(x) and of its derivatives. Thus, the extremum condition I^M- = 0, as 
a functional equation for a(x), seems untractable. The considerable complexity of the 
functional equations that determine the extremal <j(x) configurations is the source of 
all difficulties that arise in any attempt to solve the model under consideration. 

DHN found a way around this difficulty. They have discovered that the extremal 
static a(x) configurations are necessarily reflectionless [1]. Let us briefly recall the 
arguments of [1]. DHN used inverse scattering techniques [3] to express the energy 
functional £[cr] (2.19) in terms of the so-called "scattering data" associated with, e.g., 
the hamiltonian Hj, in (2.16), and thus with a(x). The scattering data associated 
with Hb are [3] the reflection amplitude r(k) of Hb (where k is the momentum of 
the scattering state), the number K of bound states in Hb and their corresponding 
energies < uj 2 < m 2 , (n = 1, • • • K), and also additional K parameters {c n }, where 
c n has to do with the normalization of the nth bound state wave function ip n . 

More precisely, The nth bound state wave function, with energy u 2 , must decay 
as ip n {x) ~ const, exp where 



If we impose that ip n {x) be normalized, this will determine the constant coefficient as 
c n . (With no loss of generality, we may take c n > 0.) 

Thus, to summarize, in the inverse scattering technique, one trades the indepen- 
dent variables a(x) for the scattering data. Then, one looks for scattering data that 
extremize S[a\. 

The key point is that £[a] depends on r(k) only through certain dispersion inte- 




(3.1) 
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grals involving log(l — |r(A;)| 2 ). Thus, the saddle point condition 




(3.2) 



is almost trivially solved by 



r(k) = 0. 



(3.3) 



Moreover, there seem to be no other solutions of (3.2). In other words, extremal 
static a(x) configurations are necessarily reflectionless. 

This restriction on cr(x) is very powerful, since once the reflection amplitude, i.e., 
the function r(k), is eliminated out of the scattering data, all that remains is a discrete 
set of 2K real parameters: 



Thus, the space of all reflectionless backgrounds o~(x) is parametrized by a discrete 
set of real parameters. We still have to extremize the energy functional £[o~] with 
respect to these parameters. 

Note, that due to the elementary fact, that the spectrum of the one dimensional 
Schrodinger operator H b cannot be degenerate, all the u; n 's must be different from 
each other. Thus, the inequalities in (3.4) are strict. (See also the remark following 



We see that the formidable problem of finding the extremal o~(x) configurations 
of the energy functional £[a] (2.19), is reduced to the simpler problem of extremizing 
an ordinary function £(cu n ,c n ) = £ [a(x;cu n ,c n )] with respect to the 2K parameters 
{c n ,u n } that determine the reflectionless background a(x). If we solve this ordinary 
extremum problem, we will be able to calculate the mass of the fermion bag. 

In Appendix B we summarized (in some detail) the necessary facts about reflec- 
tionless potentials and the diagonal resolvents associated with them. (We will need 
these resolvents below.) In particular, Eqs. (B.38) and (B.39) of Appendix B tell us 
how the 2K parameters determine o~(x) explicitly. For the particular examples with 
K — 1 and K = 2, see (B.46) and (B.64), respectively. 



< uji < u 2 < • • • < lok < m 



ci,c 2 , - • • ,ck ■ 



(3.4) 



(B.3).) 
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3.1 Extremal Static Reflect ionless Fermion Bags 

The energy functional S[a], evaluated on a reflectionless a(x) configuration, is an 
ordinary function of the parameters (3.4) which define a(x). This function (when 
defined with the renormalized £ in (2.19)) 

M(u n , c n ) = £[cr(x; u n , c n )\ (3.5) 

is the mass of the corresponding fermion bag. Thus, the extremal bags are determined 
by solving 

dM dM . , 

^— = , ^- = 0. 3.6 

duj n dc n 

Let a be any one of the 2K parameters in (3.4). Then, 

8M_~M*p., (3.7) 



da J 5a(x) da ' 

— oo 

where J^j- is evaluated at the appropriate reflectionless a(x). 

3.1.1 The Extremum Conditions on the Mass M(u n ,c n ), Its Flat Direc- 
tions, and Collective Coordinates. 

8S \cr\ 

In order to calculate around a generic a{x) we start with the elementary identity 

oo 

5Ti\og(-d 2 x + V -u 2 ) = J dx{x\ _ d2 + v _ uj2 \x) SV(x), (3.8) 

-oo x 

and apply it to (2.19), with Hb and H c defined in (2.16). For convenience, let us 
record here the corresponding potentials, 

V b = a 2 -a' and V c = a 2 + a' . (3.9) 

Thus, using (2.17) (or (A. 31)), which relate the diagonal resolvents of H b and H c , 
respectively, to the entries B(x) and C(x) in the diagonal resolvent (2.11), we arrive 
at 



. (3.10) 
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By substituting (3.10) into (3.7) we obtain 



dM 
da 



= dx 



1 da 2 



N 



2g 2 da 



duo 
2^ 



'B(x)\ dV h (x) 



'C(x)\ dV c (x) 



LO 



da 



LO 



da 



(3.11) 

As was mentioned above, the diagonal resolvent of the Dirac operator in reflec- 
tionless a(x) backgrounds was evaluated in Appendix B (see section B.2). Thus, our 
next step is to substitute B(x) and —C(x) from (B.32) in (3.11). Let us do this 
explicitly for B(x). According to (B.32) 



B(x) = (x\ 



LO 



LO 



H b - lo'< 



X 



2%/ 



m' 



LO" 



' K 

1-2E 



n=l 



LO' 



LOt 



(3.12) 



where ip n are the (normalized) bound state wave functions of H b (and we have also 
used k = i\Jm? — lo 2 ). Thus, 



OO / 



B(x)\ dV b (x) 



LO 



da 



2V' 



i 7 ( K 

^== / dx l-2£ 

171 - u JL V n=l 



K n i) 2 n \ dV b {x) 



LO' 



oot 



da 



2\f; 



m' 



LO" 



da 



n=l 



LO' 



oot 



(3.13) 



From first order perturbation theory we know that 



// \ dVb \l 



(^n\^-\fpn) 



<-» T III \ T fit r-. 

da da 
Thus, we may simplify (3.13) further and obtain 

'B(x)\ dV b (x) 1 



du 2 n 
da 



(3.14) 



dx 



LO 



da 



m' 



LO' 



J da 



n=l 



(3.15) 



Recall that H b and H c are isospectral, save, possibly, a zero energy ground state, 
which can appear in the spectrum of only one of these operators. Since the zero 
energy ground state (if it exists) is the only difference in the energy spectra of H b and 
H c , and since, in any case, it obviously does not contribute to the sum in (3.15), we 
conclude that the contribution of —C(x) to (3.11) is the same as (3.15), but with V c 
instead of V b in the first term in (3.15). 
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Thus, combining the contributions of of B(x) and —C(x) into (3.11) we obtain 



™ = (J- + iN (*!L \ ) 7 dx — -iNY I*?- 1 Kn{ ^ 

da \2g 2 ^ J 2n 2^m 2 - u 2 ) J da ^ J 2tt ^m 2 - u 2 uj 2 - uj 2 



(where in the first term we have used + V c — 2a 2 ). 

Consider now the integrals over u in (3.16). We will show in Appendix C (see 
(C.4)) that the bare gap equation may be written as 



The remaining integral over oo is UV-convergent. Thus, renormalization of the cou- 
pling constant g 2 , based on the gap equation (see (2.5)), also renders (3.16) finite. 
Next, observe from (3.17), that if a is one of the coefficients c& in (3.4), then 



identically! This does not produce any condition on the c^'s. The energy functional 
S[a] (2.19), evaluated on a reflectionless a(x; uj n , c n ), is independent of the coefficients 
Cfc that do affect the shape of a(x). The c^'s are thus flat directions of S[a] in the 
space of all reflectionless o~(x) configurations. In fact, we show in Appendix B, that 
the Cfc's (or more precisely, their logarithms) are collective translational coordinates 
of the fermion bag a(x) (see e.g., (B.46) and (B.64)), which is a familiar concept in 
soliton and instanton physics. One of these coordinates, corresponds, of course, to 
global translations of the bag as a whole. 

The full effective action functional S e ff[a(x, t)] contains, in principle, all informa- 
tion on dynamics of space-time dependent a(x, t) bags, or solitons. However, it is 
a complicated and almost intractable object. Some progress may be achieved, per- 
haps, by trying to extend the periodic time- dependent solution of (2.3) which was 



(3.16) 




Thus, the first term in (3.16) vanishes, and 




(3.17) 
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guessed by DHN in [1]. As another small step towards understanding time dependent 
soliton dynamics from S e ff, one might consider elevating the c^'s, which determine 
the shape of static reflectionless bags, into slowly varying functions of time, and thus 
study soliton dynamics in the framework of an adiabatic approximation, in the spirit 



Let us return to (3.17), and consider the remaining possibility a = u n . In this 
case we obtain 



which has to vanish for extremal bags, and thus produces an equation to determine 
uj n . Note that (3.19) involves only one particular uj n . Thus, the extremal u; n 's are 
determined independently of each other, except for the constraint that they cannot 
coincide. 

We observe from (3.19) that ui = is a possible solution of 6 J^p = 0. Thus, there 
are extremal bags for which H b (or H c ) has a zero energy ground state. As explained 
in Appendix A (see subsection A. 1.1), a zero mode u\ = occurs as a solution of 
(3.19) if and only if cr(x) has non-trivial topology. All other solutions of = 
must be strictly positive. 

In our calculations so far we used DHN's result that extremal static a(x) con- 
figurations are necessarily reflectionless, and looked for the extremal reflectionless 
configurations. As a consistency check of our calculations, we verify in Appendix 
C that vanishing of the right hand side of (3.19) guarantees that the corresponding 
reflectionless a(x) is indeed extremal among all possible static configurations. 

3.1.2 Solution of the Extremum Conditions and Computation of the Soli- 
ton's Mass. 

In order to determine the extremal cr(x) configurations, we have to determine the 

UnS from (3.19). With no loss of generality we will assume henceforth that u n > 0. 

6 According to calculations in the next subsection, the integral J 2 ^" a ^== ^ 2 (along the appro- 
priate contour) is finite. Thus, the right-hand side of (3.19) vanishes at u>\ = 0. 



of [28, 29]. 




(3.19) 
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In order to determine the non-vanishing cj„'s, we have to evaluate the (UV finite) 
integral 



To this end we have to choose the proper contour C, and thus we have to invoke 
our understanding of the physics of fermions: The Dirac equation (i$ — a(x))^ = 
(Eq. (1.6)) in the background a(x) under consideration has a pair of charge-conjugate 
bound states at ±u n . These are the simple poles in (3.20). The continuum states 
appear as the two cuts along the real axis with branch points at ±m (see Figure 1). 

The bound states at ±oo n are to be considered together due to the charge conju- 
gation invariance of the GN model, as we discuss in Appendix D (see in particular 
subsection D.2.2). Due to Pauli's principle, we can populate each of the bound states 
±u n with up to N (non-interacting) fermions. In such a typical multiparticle state, 
the negative frequency state is populated by N — h n fermions and the positive fre- 
quency state contains p n fermions. In the parlance of Dirac's hole theory, we have 
thus created a many fermion state, with p n particles and h n holes occupying the pair 
of charge-conjugate bound states at energies ±u n . Let us name this many body state 
a (p n , hn) configuration. 

Since in this paper we are interested only in solitons in their ground state, we must 
take all states in the Dirac sea completely filled, and no positive energy scattering 
states. 

Mathematically, we thus have to let C enclose the cut on the negative u axis N 
times, and then go iV — h n times around the pole at — u n , and p n times around the 
pole at u n , as shown in Figure 1. 




(3.20) 
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m 



Fig. 1: The contour C in the complex u> plane in Eq. (3.20). The continuum states appear as 
the two cuts along the real axis with branch points at ±m, and the bound states are the poles at 
±u) n . The contour wraps N times around the cut on the negative u) axis, and then go N — h n times 
around the pole at — u n , and p n times around the pole at w n . 



Thus, the integral (3.20) depends also on h n and p n . In this way, we obtain 7 



(with n n = ym 2 — u%). The first term in (3.21) is the contribution of the part of C 
wrapped around the negative cut, and the other two terms come, respectively, from 
the poles at u) n and — u n . The function I(cu n ,p n ,h n ) depends on h n and p n only 
through their sum 



the total number of particles and holes in the pair of charge conjugate bound states 
at ±uj n . Pauli's Principle allows all < v n < 2N. However, it turns out that the 
allowed range for u n in extremal configurations is restricted to only half that range: 
Substituting (3.21) into (3.19) we obtain 




(3.21) 




V n =Vn + h 



(3.22) 




(3.23) 



Thus, the extremal value a;* satisfies 




(3.24) 



rendering 




(3.25) 



This expression, for h n = 0, is quoted, e.g., in [1, 16]. 
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Thus, since u;* > 0, we must have 8 v n < N . 

However, the end-points of this range cannot be realized physically. Concerning 
u n = 0, note that if the pair of energy levels corresponding to (3.25) are empty, then 
uj* n = m plunges into the continuum and ceases to be a bound state. Thus, the self- 
consistent, extremal energy levels (3.25) must trap fermions in order to exist. This 
fact will be related to soliton stability in Section 4 (see also the discussion around 
(3.32) in this section). At the other end-point v n = N, we obtain a>* = 0. As we have 
already mentioned (recall the discussion in subsection A. 1.1 in Appendix A, which 
we have already referred to following (2.16) and (3.19)), = occurs if and only 
if a(x) has non-trivial topology, and does not depend on dynamical details such as 
the value of v n . Thus, v n = N cannot occur in a topologically trivial soliton. What 
is perhaps less obvious, is that we cannot have v n = N in a topologically non-trivial 
soliton either. To clarify this assertion, consider an extremal topologically non-trivial 
soliton with = cj* < lo\ < ■ ■ ■ < uo* K . Then, form a sequence of such solitons in 
which z/ 2 increases. As u 2 tends to N, uj\ tends to coalesce with uj\ = 0. This cannot 
happen, because the spectrum {cu* 2 } of the one-dimensional Schrodinger operator Hf, 
(or H c ) is not degenerate. Putting all these facts together we conclude that v n is 
restricted to the range 

< v n < N . (3.26) 

We would like at this point to integrate (3.23) and find the mass M(ui, . . . ,ouk). 
To this end it is convenient to introduce [1] the angular parameter 6 n through the 
relations 

u n = mcos9 n and K n = msm6 ni (3.27) 

where positivity of u) n and n n implies < 9 n < |. In particular, note that the 
extremal value u* in (3.25) corresponds to 

0*n = ^- (3-28) 

8 This can be already seen from (3.23). From the actual computation of the contribution of the 
part of C wrapped around the negative cut, the first term in (3.21), it is clear that arctan^ > 0. 
Thus, = has a solution only in the range (3.26). 
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Recall that all the u; n 's which parametrize a given static soliton, must be different 
from each other (see the paragraph following (3.4)). Consequently, all the v n must 
be different from each other, and so must be the 0*'s. 
In terms of 9 n , we may rewrite (3.19) as 

dM 1 dM 2N6„ 



doj n msm9 n d6 n n 



(3.29) 



Finally, we can integrate (3.29) back, to obtain M(9 1 , ■ ■ -6 K ) explicitly. Thus, 



K 



M(8 ir --,8 K )=m 



71 = 1 



2N6A 2iV 
v n cos v n H sin u % 

It / 7T 



+ const. (3.30) 



Note that the mass (3.30) depends only on the # n 's which correspond to strictly 
positive u>nS. The integration constant in (3.30) depends on the topology of a(x), 
and will be determined below according to physical considerations. 

Eq. (3.30) is the mass of a fermion bag with an arbitrary reflect ionless profile a(x). 
We are interested only in extremal configurations. Substituting (3.28) in (3.30), we 
obtain the mass of the extremal bag as 

M(9t, ■ ■ ■ , P K ) = ^ E sin (^) + const. (3.31) 

oj n >0 

The extremal value (3.25) of the parameter u n is determined by the total number 
u n of particles and holes trapped in the bound states of the Dirac equation at ±u n , 
and not by the numbers of trapped particles and holes separately. As we explain 
in Appendix D (see subsection D.2.2), this fact is a manifestation of the underlying 
0(2N) symmetry, which treats particles and holes symmetrically. It indicates that 
this pair of bound states gives rise to an 0(2N) antisymmetric tensor multiplet of 
rank v n of soliton states. The states in this multiplet correspond to all the possible 
(p n , h n ) configurations, with v n = p n + h n fixed, subjected to Pauli's principle. There 
are T l ' / k =o C N C N~ k = C 2N = (2N)\/uJ(2N - v n )\ such states, precisely the dimension 
of an 0(2N) antisymmetric tensor of rank v n . Due to the restriction (3.26) only 
tensors of ranks < v n < N are allowed. 
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The antisymmetric tensors of higher ranks (in the range N < v n < 2N), each of 
which is dual to an antisymmetric tensor in the allowed range (3.26), are not realized 
as extremal fermion bags (their ranks cannot be interpreted as the total number of 
particles and holes trapped by the bag). Note, however, as a mathematical fact, that 
the mass formulas (3.30) and (3.31) are invariant against the duality transformation 
of tensorial ranks v n — > 2iV — u n , performed simultaneously with Q* n — > 7r — 0*. 
This is equivalent to u* n — > — u* n . Thus, a (p n , h n ) configuration belonging to an 
antisymmetric tensor multiplet of rank < u n < N is mapped in this way onto 
an (N — h n , N — p n ) configuration belonging to the dual antisymmetric tensor of 
rank 2iV — v n . Evidently, any two dual antisymmetric tensor multiplets have the 
same dimension, and make identical contributions to the mass formulas. Thus, these 
duality transformations may be used to extend the domain of definition of (3.30) and 
(3.31) to all 0(2N) antisymmetric tensors < v n < 2N . 

The bound state at uj\ = 0, which arises if and only if a(x) has non-trivial topology, 
behaves differently from the previous case. Due to its topological origin, this mode 
is stable and cannot flow away from zero as the other parameters which define a(x) 
vary. In particular, it remains unchanged as the number of trapped fermions varies, in 
contrast with (3.25). Since this bound state is unpaired, it can trap up to N fermions. 
Thus, it gives rise to a 2^ dimensional multiplet of soliton states, which, as we show 
in subsection D.2.1, is the spinorial representation of 0(2N) [30]. 

Let us determine now the integration constant in (3.31). As was mentioned earlier, 
this integration constant depends on the topological charge carried by the correspond- 
ing extremal cr(x) configuration. We have already mentioned in subsection 2.1, that 
topologically trivial solitons are stable because of the binding energy released by the 
trapped fermions, and cannot form without such binding. Thus, the binding energy 

B(9l, ...,e* K )=mJ2"n- M(0l ■■■,9* K ) (3.32) 

n=l 

of a topologically trivial soliton, which measures the stability of the soliton, should 
tend to zero as the number of trapped fermions tends to zero. In other words, if each 
of the pairs of levels ±u n traps a small number u n « N of fermions, we expect that 
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M(6>i, • • • , 9* K ) ~ rnY,n=i v n- Imposing this physical constraint on (3.31) determines 
the integration constant (in the topologically trivial case) to be zero. This cannot 
be true in the topologically non-trivial sector, since the the integration constant in 
(3.31) is evidently the minimal mass in that sector. The latter is the mass of the 
topologically non-trivial soliton which has a single bound state at lo\ — 0, namely the 
CCGZ kink (or antikink). Its mass is 

Nm 

M kink = , (3.33) 

71 

as quoted in [8] (see also Eq. (3.30) in [1]), and is independent of the number n 
of fermions trapped in the bound state at ui = 0, as we mentioned in the previous 
paragraph. In Appendix E we will present our own derivation of (3.33), based on 
(2.19) and on a simple dimensional argument. In addition, we will provide below (see 
the passage following (3.36)) a simple physical argument for the correctness of (3.33). 
Thus, gathering all these facts together, we conclude that 

MK ... lUK ;q) = — £ sin ^ j + — * M>1 , (3.34) 

uj n >0 

where q is the topological charge (2.9) carried by the soliton (and all the ranks v n are 
different from each other). 

From (3.33) and (3.34) we see that the kink, as well as other extremal solitons 
(whose ranks v n are a finite fraction of N), are heavy objects. Their masses are of 
0(N). In fact, this should be expected, based on our experience with soliton physics 
in weakly interacting field theories: in this paper we study the GN model (1.1) in the 
large N limit, in which N ~ g~ 2 , and thus all soliton masses are of O (g~ 2 )- 

3.1.3 Summary: Masses, Profiles and 0(2N) Quantum Numbers of Ex- 
tremal Fermion Bags 

We complete the task of determining the extremal static soliton configurations by 
combining the results of the previous discussion and the results of the analysis made 
in sections B.2, B.3 and D.2 in the appendices. We summarize our findings as follows: 
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(a) topologically trivial fermion bags, q = 
An extremal static, topologically trivial soliton, bearing K pairs of bound states, 
has its K cj„'s in (3.4) determined by (3.25), with the K moduli c n left as arbitrary 
parameters. Its mass is given by (3.34) with q = 0. Its profile <j(x) is computed in 
Appendix B (see section B.2) and is given by (B.38). According to Appendix D (see 
subsection D.2.2) we should refer to this object more precisely as an 0(2N) multiplet 
of solitonic states. This multiplet transforms under 0(2N) as the direct product of 
the K individual antisymmetric tensors associated with the K a;„'s. 

In particular, the case K — 1 with U\ = mcos (f^) > 0, in which the pair of 
bound states at ±lui trap v\ — v fermions, is the so-called DHN soliton. Its mass is 

"M-^Q (3 - 35) 

and its profile is given by (B.46): 

1 , / m + k 



1 fm + K\ 

KX+- log 

2 \m — k J 



(3.36) 



a(x) = m + /ctanh(/cx) — Ktanh 

(where we have set <r(oo) = m and x = in (B.46)), with k = msin (^)- These 
formulas agree with the results quoted in Eqs. (3.27) and (3.28) of [1]. The profile 
(3.36) is that of a bound state of a kink and an antikink, with interkink distance 
2^ log (^z^)- This interkink distance diverges in the limit k — > m, i.e., when v — > N, 
the inaccessible end-point of (3.26). Heuristically, we can understand this diver- 
gence as a necessary compromise between dynamics and topology: As v tends to N, 
uj = m cos (fi^) tends to zero, which is not in the spectrum of the corresponding 
Dirac operator i$ — a(x), since a(x) is topologically trivial. The best our soliton can 
do is to split into infinitely separated kink and antikink. Probing any finite neigh- 
borhood around either of these defects may lead us to the illusion that we are in the 
topologically non-trivial sector. 

One might perhaps interpret this as a situation in which the infinitely separated 
kink and antikink are bound at threshold. (This is subjected to the plausible assump- 
tion that there are no long-range intersoliton forces in this model, since the mass gap 
is finite, and there are no gauge interactions.) Thus, (3.35) must tend, in the limit 
v — > N, to 2M kink , from which we see that M kink = ^p, in accordance with (3.33). 
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(b) topologically nontrivial fermion bags, q — ±1 
An extremal static, topologically non-trivial a(x) configuration, with 2K parameters 
(3.4), has an unpaired bound state at lo\ = 0, and additional K — 1 pairs of bound 
states at ±a> n , which are determined by (3.25). Its mass is given by (3.34) with \q\ = 1, 
and its profile is determined as plus or minus the profile given by (B.39) (depending 
on the topological charge q it carries). As explained in detail in Appendix D (see 
subsections D.2.1 and D.2.2), it gives rise to a multiplet of solitonic states which 
transform under 0{2N) as the direct product of one factor of the 2 N dimensional 
(reducible) spinor representation of 0(2N) (associated with the bound state at u\ = 
0), and the K — 1 individual antisymmetric tensors associated with the remaining 
non-vanishing o; n 's. In particular, the case K = 1 corresponds to the CCGZ kink 
(or antikink) with mass (3.33) and profile a(x) = ±mtanh (mx) (see (B.49)). These 
kinks are thus pure 0{2N) isospinors [30]. 

3.2 Fermion Number Content of the Solitonic Multiplets. 

Expressions for the expectation values of fermion bilinear operators in a static a(x) 
background are derived in Appendix A (see Section A. 3). They are given as dispersion 
integrals over linear combinations of the entries of the diagonal resolvent (2.11). In 
particular, according to (A. 43), (A. 44) and (A. 48), the regularized fermion density in 
a given background a(x) is 

0'°(*)>re g = iNf^ [(B(x) - B VAC ) - (C(x) - C VAC )\ , (3.37) 

where we have subtracted the UV-divergent vacuum contributions. 

For evaluating {j°(x)) veg in reflect ionless a(x) backgrounds, we will substitute, 
B(x) and — C(x) from (B.32) (or equivalently, the expression (3.12) for B(x), and its 
analog for —C{x)) in (3.37). For extremal reflectionless a(x) backgrounds, the focus 
of our interest in this paper, the non- vanishing w n 's are determined by (3.25). Thus, 
for each u n > 0, we will assume that the fermions which occupy the pair of energy 
levels ±u n form a (p n , h n ) configuration, with p n + h n = u n , in accordance with (3.22). 
If u>i — (i.e., if a(x) has non-trivial topology), we assume that this level is occupied 
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by < no < N fermions. 

Let us evaluate the contribution of the S-term in (3.37) explicitly. From (3.12) 
and (A. 44) we obtain this contribution as 

, K f B \ ■ < 3 - 38 > 

^ Jc„ 2m ^m 2 - u 2 to 2 - u 2 
where ipBn is the nth bound state wave function of Ht>, and C n is the contour in Fig. 
1 (with the obvious exception, that if uj\ = 0, then the contour C\ wraps n times 
around the pole at u — 0). For each uj n > 0, the contribution to (3.38) coming 
from the integral around the left cut in Fig. 1 is — y^s n , the contribution from the 
pole at — uj n is N ~ 2 hn ^> Bn , and that from the pole at u n is ^- tfj Bn . Adding these three 
terms, we obtain the total contribution of the pair of levels at ±u n to (3.38) simply as 
Pn ~ hn ip Bn . If <j(x) has non-trivial topology, and carries topological charge q — 1, then 
according to the discussion following Eqs. (A. 18), (A. 19), only has a normalizable 
zero mode ipso , and we must set = in (3.38). In this case, the contribution to 
(3.38) is (— y + ^o) ^boj where the first term comes from the integral around the 
left cut in Fig. 1 (as before), and the second term comes from the pole at uj = 0. 
Combining all these results together, we obtain that the B-term in (3.37) yields 

(-y + no) i? m S qA + £ . (3.39) 

Wn>0 

Similarly, since — C(x) is given by an expression analogous to (B.32), the C-term in 
(3.37) yields 

-y + no) Vi + E , (3-40) 



n = l 
w n >0 



where is the nth bound state wave function of H c . (Recall from the discussion 
following Eqs. (A. 18), (A. 19) that if the topological charge q — — 1, only H c has 
a normalizable zero mode ipco-) Finally, combining (3.39) and (3.40) together, we 
obtain 

0'°^) rcg =(-f + n ) (^(^ 9ll +4(^ 5 ,-l)+E (Pn-K) 

(3.41) 
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Integrating over x, we obtain the expectation value of the total fermion number Nj 
in the background of the extremal fermion bag, simply as 

(N f ) = (n - y J 6 m + E (p n - K) . (3.42) 

u) n >0 

The terms in (3.41), associated with the positive u; n 's, have simple physical interpre- 
tation: 

N^li =Vn~K = v n - 2h n , (3.43) 

the number of particles minus the number of holes, is the valence fermion number 
of the (p n , h n ) configuration, which occupies the pair of bound states at ±u n , and 
\ i^Bnix) + ^Cnix)) * s the probability density to find any of these fermions in a small 
neighborhood of the point x. The (p n ,h n ) configuration is a member of an 0(2N) 
antisymmetric tensor multiplet of rank v n < N (recall (3.26)). Thus, as we scan 
through all states in this multiplet, we see that N^ al has a symmetric spectrum 

-v n < Ng al < v n , (3.44) 

in accordance with charge conjugation invariance. 

Interpretation of the first term in (3.41), associated with the zero mode, is more 
delicate, and exhibits an interesting physical phenomenon. As in the previous case, 
(ipBo( x ) $q,i + "tPcoix) &q ,-i) is the probability density to find a fermion, trapped in the 
zero mode bound state, in a small neighborhood of the point x. However, unlike the 
previous case, the coefficient of this probability density is the valence number hq of 
fermions trapped in the zero mode, shifted by the contribution — y of the vacuum 
(the filled Dirac sea). Thus, the fermion number associated with the zero mode is 

Nf = no - y ■ (3.45) 

This quantity coincides with the fermion number operator j /\rj spmor ) (D.16) which we 
constructed explicitly for the 0(2N) spinor representation in Appendix D. 

Note that for TV odd, ivj 0) is not an integer! Eq.(3.45) thus exhibits fractional 
fermion number [31], a phenomenon which occurs because of the non-trivial topology 
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of the background interacting with the fermions, and is independent of other details 
of the background (such as a(x) being reflectionless). Furthermore, it is valid for 
any value of N, and has nothing to do with the large N limit. Charge conjugation 
invariance of the GN model restricts 9 the fractional part of to be 10 either 
or |. Indeed, it is as defined in (3.45) which acquires a symmetric spectrum 
~Y < < y, as required by charge conjugation invariance, and not the valence 
piece Uq alone. For more details on fractionalization of fermion number in quantum 
field theory, see [33]. 

9 The difference of any two eigenvalues of must be an integer, of course, and charge conju- 
gation invariance implies that if n is an eigenvalue, so is —n. 

10 In other soliton bearing quantum field theoretic models, fermion number induced by solitons 
may acquire other rational, or even irrational values [32] . 
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4 Investigating Stability of Extremal Static Fermion 
Bags 



The extremal static soliton multiplets which we encountered in the previous section, 
correspond, in the limit N — > oo, to exact eigenstates of the hamiltonian of the GN 
model. However, at large but finite N, we expect some of these states to become 
unstable and thus to acquire small widths, similarly to the behavior of baryons in 
QCD with a large number of colors [34]. The latter are also solitonic objects and are 
analogous to the "multi-quark" bound states of the GN model. 

Furthermore, we can imagine perturbing the GN action (1.1) by a small perturba- 
tion, which is a singlet under all the discrete and continuous symmetries of the model 
(e.g., by adding to (1.1) a term e J d 2 xa 2n ), and ask which of the extremal fermion 
bags of the previous section are stable against such perturbations. 

Under these circumstances, all possible decay channels of a given soliton multiplet 
must conserve, in addition to energy and momentum, 0(2N) quantum numbers and 
topological charge. 

It turns out that non-trivial results concerning stability may be established with- 
out getting into all the details of decomposing 0(2N) representations, by imposing a 
simple necessary condition on the spectrum of the fermion number operator Nf in the 
multiplets involved in a given decay channel. As we have learned so far, a given static 
soliton multiplet is a direct product of 0(2N) antisymmetric tensors and, for topolog- 
ically non-trivial solitons, a factor of the spinor representation. The decay products of 
this soliton also correspond to a direct product of antisymmetric tensors and spinors. 
We have also learned that the spectrum of Nf in the antisymmetric tensor and spinor 
representations is symmetric, namely, —Nj 1 ^ < Nf < NJ iax . When we compose 
two such representations D±, D 2 , the spectrum of Nf in the composite representation 
Di®D 2 , will obviously have the range -Nf ax {D 1 ) - Nf ax (D 2 ) < N f {D 1 <g> D 2 ) < 
Nf nax (Di)+Nf iax (D 2 ). In particular, each of the possible (i.e., integer or half-integer) 
eigenvalues in this range, will appear in at least one irreducible representation in the 
decomposition of D\ ® D 2 . More generally, the spectrum of Nf(D 1 ® D 2 - ■ • <g) D L ) 
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will have the range |N/(£>i ® • • • <8> D L )\ < Nf ax {D 1 ) + • ■■Nf ax (D L ). 

Consider now a decay process, in which a parent static soliton, which belongs to 
a (possibly reducible) representation -D parent , decays into a bunch of other solitons, 
such that the collection of all irreducible representations associated with the decay 
products is {Di, . . . , D^} (in which a given irreducible representation may occur more 
than once). By 0(2N) symmetry, the representation -D par ent must occur in the de- 
composition of Dx <S> D 2 • • ■ <S> D L . Thus, according to the discussion in the previous 
paragraph, if this decay process is allowed, we must have 

Nf aX (^parent) < N^^D,) + ■■■ + Nf aX {D L ) , (4.1) 

which is the necessary condition for 0(2N) symmetry we sought for. (Obviously, 
similar necessary conditions arise for the other N — l components of the highest weight 
vectors of the representations involved.) The decay process under consideration must 
respect energy conservation, i.e., it must be exothermic. Thus, we supplement (4.1) 
by the requirement 

M parcnt > Yl M * ( 4 - 2 ) 

products 

on the masses Mj of the particles involved. 

For each of the static soliton multiplets discussed in the previous section, we 
will scan through all decay channels and check which of these decay channels are 
necessarily closed, simply by requiring that the two conditions (4.1) and (4.2) be 
mutually contradictory. 

The oscillating time dependent DHN solitons [1] may appear as decay products in 

some of the channels. However, for a given assignment of 0(2N) quantum numbers 

of the final state, they will occur at a higher mass than the final state containing 

only static solitons 11 . Thus, if a decay channel into a final state containing only 

static solitons is necessarily closed (in the sense that (4.1) and (4.2) are mutually 

n In [1], DHN showed that the oscillating solitons occur in 0(2N) hypermultiplcts of antisymmetric 
tensors. Such a hypermultiplet is determined by a principal quantum number v p (0 < v p < N) , which 
determines the common mass sin (-^r) of all the states in it. For v v even, the hypermultiplet 
contains antisymmetric tensors of ranks 0, 2, . . . v p (each with multiplicity 1). For v v odd, it contains 
tensors of ranks 1, 3, . . . , v p . 



32 



contradictory), then all decay channels into final states with the same 0(2N) quantum 
numbers, which contain time dependent solitons, are necessarily closed as well. Thus, 
it is enough to scan only through decay channels into final states made purly of static 
solitons 12 . 

For all the DHN multiplets (K = 1, q = 0) and for the CCGZ kink multiplet 
(K — 1, \q\ — 1), we will find in this way that all decay channels are necessarily closed, 
thus establishing their stability. That these are stable multiplets is well known, of 
course [1,8]: they are the lightest solitons, given their 0(2N) and topological quantum 
numbers. Our new result is the marginal stability of the heavier, topologically non- 
trivial solitons, with a single pair of non-zero bound states (K — 2, \ q\ — 1) [23]. 
For all other static soliton multiplets, there are always decay channels in which (4.1) 
and (4.2) are mutually consistent. In such cases, more detailed analysis is required 
to establish stability or instability of the soliton under consideration. Nevertheless, 
given that the heavier topologically non-trivial soliton (K — 2, \q\ — 1) is already at 
the threshold of stability, we conjecture that all these other multiplets are unstable. 

4.1 Investigating Stability of Topologically Trivial Solitons 

Consider the decaying parent soliton to be a topologically trivial static soliton with K 
pairs of bound states at non-zero energies, corresponding to the direct product of K 
antisymmetric tensor representations of ranks i>i, . . . , vk- The corresponding angular 
variables Oi, . . . , 9^ are determined by (3.28), namely, O n = Thus, according 
to (3.34) and (3.44), the mass of this soliton is M(z> 1? . . . , v K ) = Y^n=i s i n @n, and 
the maximal fermion number eigenvalue occurring in this representation is N™^ (D pa , rent ) = 

J2 n =l ^n- 

Following the strategy which we laid above, we shall now scan through all imagin- 
able decay channels of this parent soliton (into final states of purely static solitons), 
and identify those channels which are necessarily closed. 

12 This statement is surely true if DHN's time dependent solitons exhaust all the stable time 
dependent configurations of the model. However, even if there arc time dependent solitons in the 
GN model beyond those discovered by DHN, it is most likely that they will also be heavier than the 
corresponding static states with the same 0{2N) quantum numbers. 
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We first scan through all decay channels into topologically trivial solitons. Thus, 
assume that the parent soliton under consideration decays into a configuration of 
lighter topologically trivial solitons, with quantum numbers of the direct product 
of L antisymmetric tensor representations i/i, . . . ,ul (and corresponding angular pa- 
rameters 0i, ... , 0l). The way in which these L multiplets are arranged into extremal 
fermion bags is of no consequence to our discussion. Thus, we discuss all decay 
channels consistent with these quantum numbers in one sweep. 

The necessary conditions (4.1) and (4.2) imply 

K L K L 

E e "<E^ E sin0 n>E sin ^> ( 4 - 3 ) 

n=l i=l n=l i=l 

where all < 6 n , 0i < |. The two pairs of boundary hypersurfaces in (4.3) are 

Ei, Si: e 1 + ... + e L = e 1 + ... + e K 

E 2 ,E 2 : sin 1 H h sin 0^ = sin @i H hsin©^, (4.4) 

where Ei ;2 are hypersurfaces in #-space and Ei i2 are the corresponding hypersurfaces 
in O-space. 

The two necessary conditions in (4.3) will contradict each other, thereby protecting 
the parent soliton against decaying through the channel under consideration, if and 
only if the hypersurface S 2 lies (in the positive orthant) between the origin and the 
hyperplane Si. When will this happen? 

To answer this question it is useful to split 

sinOi + h sin ©^ = s + sin f3 , where (4.5) 



s = integral part of [sin 6i + • • • + sin Q K ] and 
sin p = fractional part of [sin ©i + • • • + sin ©x] • (4.6) 

Thus, the integer < s < K, and < (3 < f . 

The variables 0i are restricted to the L dimensional hypercube [0, |] L in the 
positive orthant. The hypersurface S 2 intersects the boundaries of this cube along 
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(curved) subsimplexes whose vertices are 

0?° = \ + fe 2 + • • • + S iis ) + (36 iis+1 , (i = 1, . . . , L) (4.7) 

with all possible choices of s + 1 coordinates i±, . . . ,i s +i out of L (as the reader can 
convince herself or himself by working out explicitly a few low dimensional examples). 
All these vertices lie on the hypersurface 

S s/3 : 9 1 + --- + e L = /3+^s. (4.8) 

The curved hypersurface E 2 is dented towards the origin. In the positive orthant it 
lies between the hypersurface E s/ 3 and the origin. Thus, the two necessary conditions 
in (4.3) will contradict each other as required, if and only if E s/3 lies between Ei and 
the origin, namely, when the parameters of the parent soliton are restricted according 
to 

&i + --- + &k>P+^s. (4.9) 

At the same time, these parameters are subjected to (4.5), which picks up a particular 
hypersurface E 2 in the space of the 6's. Thus, the desired contradiction of the two 
necessary conditions in (4.3) occurs for parent solitons which correspond to points on 
the intersection of E 2 and (4.9). Let us determine this intersection: 

The hypersurface E 2 intersects the boundaries of the K dimensional hypercube 
[0, |] K (the range of the 6 n 's) along (curved) subsimplexes whose vertices are 

©n } = |(<W + Snn 2 + ■ ■ ■ + Snn s ) + PSnn a+1 , ijl = 1, . . . , K) (4.10) 

with all possible choices of s + 1 coordinates ii, . . . ,i s +i out of K (similarly to the 
intersection of E 2 and the boundary of the L dimensional cube [0, |] L ). All these 
vertices lie on the hypersurface 

s s/3 : e 1 + --- + e K = p + ^s, (4.ii) 

which is, of course, the boundary of (4.9). Furthermore, E 2 is dented towards the 
origin and lies (in the positive orthant) between the origin and E S/ g. Thus, we conclude 
that the desired intersection is precisely the set of all vertices defined by (4.10). 
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Each such vertex represents a soliton which cannot decay through the channel 
under consideration, and is thus potentially stable. More precisely all these vertices 
correspond to the same soliton, since the coordinates of these vertices are just permu- 
tations of each other, and thus all of them correspond to the same set of parameters, 
in which 

7T 

s of the 6's are degenerate and equal to — 
one of the O's is equal to (5 , and 

the remaining K — (s + 1) O's are null. (4-12) 

In Appendix F we provide an alternative proof of (4.12). 

Does such a soliton exist? To answer this question let us recall a few basic facts: 
The parent soliton under discussion is topologically trivial. As such, it must bind 
fermions to be stabilized, and none of its bound state energies may vanish. Thus, all 
the ranks occurring in it must satisfy < i> n < N, in accordance with (3.26). Finally, 
recall that all the o; n 's which parametrize a given static soliton, must be different 
from each other (see the paragraph following (3.28)). Gathering all these facts we 
deduce from (3.27) and (3.28) that all the 9 n 's must be different from each other, 
and furthermore, that each 6 n ^ 0, |. 

Thus, the only physically realizable parent solitons, which are necessarily stable 
against the decay channel in question, correspond to s = and K — s + 1 — 1. These 
are, of course, the static DHN solitons with profile (3.36) and mass (3.35). 

A corollary of our analysis so far is that DHN solitons are stable against decaying 
into Lf free fermions or antifermions 13 (i.e., Lf fundamental 0(2N) representations) 
plus L — Lf solitons corresponding to higher antisymmetric tensor representations. 
In particular, it is stable against complete evaporation into free fermions. This fact 
manifests itself in the binding energy function (3.32), which for a DHN soliton of rank 

TT 

B(e) = ^—(e-sme). (4.13) 

7T 



13 Strictly speaking, this argument is valid only for values of Lf which arc a finite fraction of N, 
since our mass formula (3.34) is the leading order in the 1/N expansion, while removing a finite 
number of particles from the parent soliton is a perturbation of the order 1/N relative to (3.34). 
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B(Q) is positive and increases monotonically in the physical range < < f , and 
so does the binding energy per fermion, — m (l — • This is in contrast 
with binding in nuclei, where the binding energy per nucleon saturates. The DHN 
soliton becomes ever more stable as it traps more fermions (up to the maximal number 
v = N). This is the so-called "mattress effect", familiar from the physics of fermion 
bags. 

All the DHN solitons (K = 1, q = 0) have masses below the kink-antikink thresh- 
old Thus they cannot decay into pairs of topologically non-trivial solitons. The 
topologically trivial solitons with K > 2, on the other hand, may have masses above 
the kink-antikink threshold, but there is not much interest in looking for such heavier 
solitons which are stable against decaying into pairs of topological defects, since we 
have already established that they are not necessarily protected against decaying into 
topologically trivial solitons 14 . 

To summarize, we have thus covered all possible decay channels of a given static 
topologically trivial soliton into a configuration of other solitons (static or time de- 
pendent). We have found that the only static solitons, which are necessarily stable 
against decaying into any of these channels (in the sense that (4.1) and (4.2) will al- 
ways contradict each other), are just the DHN solitons. This is of course, well known. 
It follows almost trivially from the fact that if < 9 n < |, Vn and < J2n@n < f, 
then sin (£„ 9 n ) < EnSin6»„. 

4.2 Investigating Stability of Topologically Non- Trivial Soli- 
tons 

The CCGZ kink and antikink are the lightest states with topological charge q = ±1. 

Thus they are stable. Are there any other stable topologically non-trivial solitons? 

To answer this question we will apply our strategy to specify all topologically 

non-trivial solitons which are necessarily stable. Thus, consider such a generic parent 

soliton. It has a zero energy bound state iQ\ — with its corresponding spinorial 

14 For what it worths, it can be shown, by going through the same analysis as in the previous 
case, that there are no decay channels into pairs of topological solitons, against which any of the 
K > 2, q = solitons (with mass above the kink-antikink threshold) is necessarily protected. 
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representation factor, and additional K — 1 pairs of bound states at non-zero u; n 's, 
with the direct product of their K — 1 antisymmetric tensor representations of ranks 
z>!, . . . , z>^_i. The corresponding angular variables 61, . . . , Qr-i are determined by 
(3.28), namely, O n = 7^. Thus, according to (3.34), (3.44) and (3.45), the mass of 
this soliton is M(z>i, . . . , vk-\) = + Hn=\ sin an d the maximal fermion 
number eigenvalue occurring in this representation is N™ ax (D pa _ Te nt) = y + S^Ti 1 ^n- 

The parent soliton may decay into a final state, which will contain, in addition 
to a lighter soliton with the same topological charge, P pairs of topologically non- 
trivial solitons (with opposite topological charges), and a bunch of topologically trivial 
solitons. The final state of the decay is thus most generally a direct product of 
some L antisymmetric tensor representations Ui, . . . , vl (and corresponding angular 
parameters 6*i, . . . , 9 L ), and 2P + 1 spinorial representations. As in the previous case, 
the way in which the L antisymmetric tensors are arranged into extremal fermion 
bags is of no consequence to our discussion. 

The necessary conditions (4.1) and (4.2) now imply 

E ®n < \P + E Q i > sin 6 n > P + E sin 9 t , (4.14) 

n=l i=l n=l i=l 

(where all < n , 0j < |). These conditions are identical to (4.3), with K there 
replaced here by K — 1 and with s > P. Going through the same analysis as we did for 
the topologically trivial soliton, we conclude that the two necessary conditions (4.14) 
contradict each other only when s = (which means P = 0), and K — 1 < s + 1 = 1. 

Thus, in addition to the CCGZ kink (K — 1), we have discovered a heavier stable, 
topologically non-trivial soliton (K = 2). As in [23], we shall refer to it as the "Heavier 
Topological Soliton" (HTS), which is a bound state of a DHN soliton and a kink or 
an antikink. Its mass is 

Nm 2Nm . ( nu \ tA . 

M ffTS ,„ = — + — sm(-), (4.15) 

where v is the rank of its antisymmetric tensor factor, which is tensored with the 
2 N dimensional spinor. Thus, Mhts,u coincides with the sum of masses of a CCGZ 
kink and a DHN soliton of rank v. Drawing further the analogy between the solitons 
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discussed in this Review and baryons in QCD with large iV co / or , the HTS would 
correspond to a dibaryon. 

The profile of this HTS (for boundary conditions corresponding to q — 1) is given 
by (B.66) and (B.69), which we repeat here for convenience: 

2k 



a(x) = m + 



Y _|_ m+K e 2K(x-y ) 



Y _i_ m+K ^2m(x—xo) _i_ m+K ^2K(x-yo) 
n/ \ (m—K) 2 (m—K) 2 

— 2(m + K) 2 2 2 

I i f m+n \ e 2m(x-xo) _|_ f ™+k \ e 2K(x-y () ) i / ™+« N \ g 2m(a;-a;o)+2K(a;-2/o) 
\m—K J \m—K J \m—K J 

= —k tanh [/c(x — yo + R)] 

sinh [m(x — x ) + k(x — yo) + 2kR] + sinh [m(x — x ) — k(x — yo)] 



e kR cosh [m(x — x ) + k(x — y ) + 2kR] + e KR cosh [m(x — x ) — k(x — y )] 

(4.16) 



with R = ^log (2j±s) (Eq.(B.67)) and k = ^m 2 - cu b 2 = msin (^), as usual. The 
arbitrary real parameters a; and y are related to the coefficients Ci,c 2 in (3.4) via 
(B.65). These are the flat directions, or moduli, which we discussed in subsection 
3.1.1. We can invoke global translational invariance to set one of these collective 
coordinates, say xq, to zero. The remaining parameter yo is left arbitrary. 

In [23] we have studied the profile of (4.16) in some detail. (See Figures 1-3 
in [23].) In general, it has the shape of the profiles of a CCGZ kink and a DHN 
soliton superimposed on each other, which move relative to each other as we vary the 
remaining free modulus yo. The latter is essentially the separation of the kink and 
DHN soliton. The mass (4.15) of the HTS remains fixed throughout this variation, 
as we discussed in subsection 3.1.1. 

The fact that dM HT s,v/dyo = 0, just means that the kink and DHN soliton exert 
no force on each other, whatever their separation is. In other words - they are very 
loosely bound. Indeed, it is evident from (4.15) that the HTS is degenerate in mass 
with a pair of non-interacting CCGZ kink and a DHN soliton of rank v. If we interpret 
the latter two as the constituents of the HTS, they are bound at threshold. Thus, 
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the HTS is marginally stable [23]. 

This is a somewhat surprising result, since one would normally expect soliton- 
soliton interactions to be of the order -jj ~ N in a weakly interacting field theory. 
This is consistent, of course, with what one should expect from general counting 
rules [34], according to which the baryon-baryon interaction is of order N. Yet, due 
to dynamical reasons which ellude us at this point, the solitonic constituents of the 
HTS avoid these general considerations and do not exert force on each other. 

The limit v — > N is of some interest. Strictly speaking, there is no HTS with 
v = N, as was discussed above. However, it is possible to study HTS's with v 
arbitrarily close to N. In this case, k — > m, with R — > oo in (4.16). Thus, for 
\x\, \xq\, 1 2/o | << R, (4.16) tends in this limit to 

I _ e -2m(x-y ) _ e -2m(x-x ) 

a(x) = m- — -, r — r . (4.17) 

In the asymptotic region 1 << m\xo — y \ (<< mR) (4.17) simplifies further, and 
appears as a kink mtanh[m(x — x max )], located at x max = max{x , yo}- This clearly 
has mass M kink = ^p, but according to (4.15), M hts , v ~n should tend to = 
3Mfci n fc. The extra mass 2Mki n k corresponds, of course, to the kink- anti- kink pair 
which receded to spatial infinity. 
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Appendix A: Resolvent of the Dirac Operator With 

Static Background Fields and Super- 
symmetry 



In this Appendix we recall some useful properties of the resolvent of the Dirac operator 
D = wy° + i^ x d x — a(x) in a static background. In particular, we make connections 
with supersymmetric quantum mechanics 15 and use it to express all four entries of 
the diagonal resolvent in terms of a single function. This will be useful in evaluating 
masses of fermion bags. 

The presentation in this Appendix is a special case of the recent discussion in 
[20], which analyzes the more generic Dirac operator D = i$ — (cr(x) + iir(x)^ 5 ). 
The spectral theory of the latter operator gives rise to an interesting generalization 
of supersymmetric quantum mechanic which was discussed in [20]. We also mention 
that supersymmetric quantum mechanics was instrumental in the analysis of fermion 
bags in [16] . 

As was discussed in Section 2, our task is to invert the Dirac operator (2.10), 
D = wy° + irfdx — <j(x), in a static background a(x) subjected to the boundary 
conditions (2.8). We emphasize that inverting (2.10) has nothing to do with the large 
N approximation, and consequently our results in this section are valid for any value 
of N. 

In the Majorana representation (1.2) we are using in this paper, 7 = a 2 , 7 1 — i< J 3 
and 7 5 = — 7°7 1 = <7i, the Dirac operator (2.10) is 



d x — cr —iuj \ 


I -Q 


—ioj \ 


iuj d x — cr J 


\ iuo 


-Q f ) 



where we introduced the pair of adjoint operators 

Q = a(x) +d x , Qt = a ( x ) _ dx . (A.2) 
15 For a thorough review of supersymmetric quantum mechanics see [13]. 
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Inverting (A.l) is achieved by solving 



/ -Q 


—ioj \ 


( 


\ ioo 


-Q ] ) 


\ 



—H8(x — y) 



(A.3) 



for the Green's function of (A.l) in a given background a(x). By dimensional analysis, 
we see that the quantities a, b, c and d are dimensionless. 

A.l "Supersymmetry" in a Bag Background 

The diagonal elements a(x,y), d(x,y) in (A.3) may be expressed in term of the off- 
diagonal elements as 



a(x, y) = — Q*c(x, y) , d(x, y) = -Qb(x, y) 

CO CO 

which in turn satisfy the second order partial differential equations 



(A.4) 



(QtQ-u?)^ = 

(ggt_ w 2)«£i¥) = 



-d 2 x + a(x) 2 - a\x) - oo 2 
-d 2 x + a(x) 2 + a'(x)-co 2 



b(x,y) 



c(x,y) 



CO 



= S(x-y) 



= S(x-y). 



(A.5) 



Thus, b(x,y)/oo and —c(x,y)/u are simply the Green's functions of the corre- 
sponding Schrodinger operators 

H b = Q^Q = [-dl + a{xf - a'{x)] , and H c = QQ f = [-d 2 x + a{xf + a\x) 

(A.6) 

in (A.5), namely, 

b(x, y) 6(x- y) b 2 (x)b 1 (y) + 6 (y - x) b 2 (y)b 1 (x) 



CO 



c(x,y) 



LU 



W h 



6{x-y) c 2 (x)c 1 (y) + 9 (y - x) c 2 (y)c 1 (x) 



(A.7) 



Here {bi(x) , b 2 (x)} and {ci(x) , c 2 (x)} are pairs of independent fundamental solutions 
of the two equations 

H b b(x) = u 2 b(x) and H c c(x) = uj 2 c(x) , (A. 8) 
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subjected to the boundary conditions 

&i(a0,ci(x) — A£le- ikx , b 2 (x) , c 2 (x) — A^e ite (A.9) 

with some possibly /c-dependent coefficients and with 16 

fc = Vcj 2 - m 2 , ImA;>0. (A.10) 

The purpose of introducing the (yet unspecified) coefficients A^l(k), A^ c {k) will be- 
come clear later, following Eqs. (A. 15) and (A. 16) below. The boundary conditions 
(A.9) are consistent, of course, with the asymptotic behavior (2.8) on a due to which 
both H b and H c tend to a free particle hamiltonian — <9 2 + m 2 as x — > ±oo. 
The wronskians of these pairs of solutions are 

W b (k)=b 2 (x)b' 1 (x)-b 1 (x)b' 2 (x) 



W c (k) = c^xjc'^x) - d{x)c 2 (x) . 

(A.ll) 

As is well known, W b (k) and W c (k) are independent of x. 

We comment at this point (for later reference) that in the scattering theory of the 
operators H b and H c it is more customary to consider pairs of independent funda- 
mental solutions {b\(x),b 2 (x)} and {c\(x), c 2 (x)} for which the coefficients A^\,A^ b 2 \ 
in (A.9) are k independent, because, as is well known in the literature [3], with 
such boundary conditions, the corresponding wronskians are proportional (up to a k- 
independent coefficient) to k/t(k), where t(k) is the transmission amplitude of the cor- 
responding operator H b or H c . Thus, we will refer to the wronskians of pairs of inde- 
pendent fundamental solutions with asymptotic behavior (A.9) with A^ c = A^l = 1 
as the standard wronskians W b st (k) and Wf(k). Therefore, we can express the wron- 
skians (A.ll) of pairs of fundamental solutions with asymptotic behavior (A.9) as 

W b (k) = A b 1 \k)A ( b 2 \k)W b st (k) and W c {k) = A^(k)Af\k)Wf(k) . (A. 12) 

16 We see that if Imk > 0, &i and c\ decay exponentially to the left, and b 2 and c 2 decay to the 
right. Thus, if Imfc > 0, both b(x, y) and c(x, y) decay as \x — y\ tends to infinity. 
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We remind the reader that at a bound state energy uj = ui b (or, equivalently, at the 
corresponding k b = i^Jm 2 — uj 2 ), say, of the operator H b , the transmission amplitude 
tb(k) of H b has a simple pole, t b (k) ~ l/(k — &&)■ Thus, W 6 S *(A;) ~ k/t b (k) vanishes as 

W 6 S *(A;) ~ - , A;~A; 6 (A. 13) 

near that bound state. This behavior occurs simply because at the bound state 
k = kb, bi(x) and b2(x) are both proportional to the bound state wave function 
(which decays asymptotically as exp(— \Jm? — ~u^\x\) ), and are thus linearly depen- 
dent. Thus, the zeros of W b st (k) determine the spectrum of bound states of H b . The 
non-standard wronskian W b {k) will have some extra k dependence resulting from the 
factors A^(k)A^(k) in (A. 12). Similar assertions hold, of course for H c , Wf and 
W c . 

Substituting the expressions (A. 7) for the off-diagonal entries b(x,y) and c(x,y) 
into (A. 4), we obtain the appropriate expressions for the diagonal entries a(x,y) and 
d(x,y). We do not bother to write these expressions here. It is useful however to 
note, that despite the 9 x 's in the Q operators in (A. 4), that act on the step functions 
in (A. 7), neither a(x,y) nor d(x,y) contain pieces proportional to S(x — y) . Such 
pieces cancel one another due to the symmetry of (A. 7) under x <-> y. 

The factorization (A. 6) of H b = Q^Q and H c = QQ^ into a product of the two first 
order differential operators Q and is the hallmark of supersymmetry. As is well 
known, this factorization means that H b and H c are isospectral, i.e., have the same 
eigenvalues, save possibly an unmatched zero mode which belongs to the spectrum of 
only one of the operators. 

Let us now recall how this property arises: The factorized equations 

Q ] Qb = uj 2 b and QQ f c = uj 2 c (A. 14) 

suggest a map between their solutions. Indeed, given that H b b = u 2 b, then clearly 

c(x) = -Qb(x) (A. 15) 

UJ 

is a solution of H c c = u 2 c. The factor - in (A. 15) is relevant for the case in which 
b(x) is an eigenstate of H b (with eigenvalue oj 2 ). In such a case, c(x) on the left hand 
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side of (A. 15) is an eigenstate of H c with the same eigenvalue u 2 , and furthermore, 
due to the factor ^, it also has same norm as b(x). Similarly, if H c c = u 2 c, then 

b(x) = -Q+ c(x) (A. 16) 

solves Hbb = u 2 b, and has the same norm as c(x), in case c(x) is an eigenstate of H c . 

We remark, that in the more generic case of the Dirac operator i$ — a(x) — 7 5 7r(:r) 
discussed in [20], the factors analogous to the factors ^ in (A. 15) and (A. 16) (factors 
of ^ , see [20], Eqs. (2.12)-(2.14)), arise already at the level of the differential 
equations analogous to (A. 14), and not from considerations of the norm of b(x) and 
c(x). 

One particular useful consequence of the mappings (A. 15) and (A. 16) is that given 
a pair {bi(x), b2(x)} of independent fundamental solutions of (Hb — uo 2 )b(x) = 0, 
we can obtain from it a pair {ci(x), c 2 (x)} of independent fundamental solutions 
of (H c — u 2 )c(x) = by using (A. 15), and vice versa. Therefore, with no loss of 
generality, we henceforth assume, that the two pairs of independent fundamental 
solutions {bi(x), b 2 (x)} and {ci(x), c 2 (x)}, are related by (A. 15) and (A. 16). 

The coefficients A^l(k), A^ c (k) in (A. 9) are to be adjusted according to (A. 15) 
and (A. 16), and this was the purpose of introducing them in the first place. 

It is clear (and of course well known from the literature on supersymmetric quan- 
tum mechanics), that the mappings b(x) <-> c(x) break down if either Ht> or H c has 
an eigenstate at zero energy. In case that zero mode exists, it must be the ground 
state of the corresponding operator Hb or H c , since these operators obviously cannot 
have negative eigenvalues. 

A. 1.1 Zero Modes and Topology 

Let us determine under which conditions such zero-energy eigenstates exist and what 
they look like. Assume, for example, that Hb has its ground state \bo) at zero energy. 
Our system is defined over the whole real axis, thus bound states are strictly square- 

2 

integrable. By assumption (b \Hb\b ) = (&o|Q^Q|^o) = IQI&o) = 0. Thus, Hb\b ) = 
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implies the first-order equation 



Q bo(x) = b' (x) + cr(x)b (x) = . 



(A.17) 




(A.18) 



■X 



cq{x) = — c' (x) + a(x)co(x) = namely co(x)=exp+ / a(y)dy . (A. 19) 



For a(x) subjected to the boundary conditions (2.8), the wave function (A.18) is 
square-integrable if and only if cr(oo) = — er(oo) = to, i.e., only when a(x) carries 
topological charge q — 1. Similarly, the wave function (A. 19) is square-integrable if 
and only if er(oo) = — er(oo) = —to, i.e., only when a(x) carries topological charge q = 
— 1. Thus, b and Co cannot be simultaneously square-integrable, and correspondingly, 
only one of the operators H b and H c can have a zero energy ground state in a given 
topologically non-trivial background a(x). Of course, none of these operators has a 
zero energy state if a(x) is topologically trivial, i.e., when q = 

To summarize, if a(x) carries topological charge q — 1, H c is strictly positive 
and isospectral to H b , except for the normalizable zero-mode (A.18) of H b . In the 
opposite case, where a(x) carries topological charge q = — 1, H b is strictly positive 
and isospectral to H c , except for the normalizable zero-mode (A. 19) of H c . In the 
third case, when a(x) is topologically trivial, q — 0, neither H b nor H c has a zero 
energy state, and they are strictly isospectral. 

Let us reformulate all this in terms of the Dirac equation (1.6). Assuming the 
existence of a bound state at uj = 0, ip(x,t) = u(x), we write 



This equation may be further reduced as follows: The matrix 27 1 is hermitean with 
eigenvalues ±1. (This basis independent property is manifest in the particular repre- 
sentation (1.2).) Thus, consider (A. 20) with i^/ 1 u ± (x) = ±u±(x). The corresponding 



il l d x — cr(x) u(x) = . 



(A.20) 
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solutions are u±(x) = u±(0) exp ± / a(y) dy. If q — 1, then only u_(x) = u-(0)b (x) 

o 

is normalizable. Similarly, if q = — 1, only = u+(0)co(x) is normalizable. And 

if q = 0, none of these solutions is normalizable. This is a general feature of the Dirac 
equation (1.6). Thus, as was first noticed in [31], a fermion which is Yukawa-coupled 
to a topologically non-trivial scalar field (as in (1.1)), will have a single, unpaired, 
normalizable zero energy mode. 

The "Witten index" [35] associated with the pair of isospectral operators H b and 
H c , which in this context may be defined 17 as 

Witten index = Af (H b ) - N (H C ) , (A.21) 

where Afo(H b)C ) is the number of normalizable zero modes of H btC - This index is of 
course a topological invariant of the space of hamiltonians H b and H c defined with 
configurations a(x) that satisfy the boundary conditions (2.8). Indeed, it is clear from 
the discussion above that in our context, it coincides with the topological charge (2.9) 
of the background a(x): 

Witten index = q . (A.22) 

Thus, the Witten index in our system is either ±1 or zero. 

If one of the operators H b or H c supports a normalizable zero mode, we say 
that supersymmetry is unbroken. When none of the operators has a normalizable 
zero mode, we say that supersymmetry is broken. Since for our system only one 
operator can support a normalizable zero mode, we may summarize these definitions 
by saying that for our system a null Witten index means broken supersymmetry. A 
non vanishing Witten index, which in our case can take only the values ±1, always 
means unbroken supersymmetry. 

Since the Witten index in our model coincides with the topological charge (2.9) 

q, we may rephrase the last paragraph by saying that only topologically non-trivial 

a(x) backgrounds lead to a Dirac operator with a normalizable bound state at u — 0. 

17 With this definition, we are identifying Hb and H Cl respectively, as the hamiltonians of the 
bosonic and fermionic sectors of the supersymmetric system. 
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A. 1.2 Wronskians and Isospectrality 

An interesting outcome of the isospectrality of H b and H c concerns their wronskians. 
Indeed, from the definition (A. 11), it follows for pairs of independent fundamental so- 
lutions {b 1 (x),b 2 (x)} and {ci(x), c 2 (x)} which are related through (A. 15) and (A. 16), 
that 

W c c 2 d x c x - c x d x c 2 b 2 d x b l - b x d x b 2 W b , > 

— = = C1&2 - c 2 bi = = — , lo f: . (A.23) 

LO LO LO LO 

The wronskians of pairs of independent fundamental solutions of H b and H c are equal 
for all lo 7^ ! Eq. (A.23) will play an important role in establishing useful properties 
of the diagonal resolvent of the Dirac operator in the next subsection. 
With no loss of generality, we may choose 

4 1 ) = A ( b 2) = 1 (A.24) 

in (A.9). The coefficients A^\A^ are then determined by (A.15): 

erf— oo) — ik 



LO 



A? = (A . 25) 

LO 

Thus, from (A.12), (A.24), (A.23) and (A.25) we have 

w b (k) = w-(k) = K-°°)-yoo) + »fc) w « = Wc{k) 



LO" 

,2 — ^2 , u2 



Using lo = m 2 + k 2 , the definition of topological charge q = (er(oo) — er(— oo))/2m 
(2.9), and the boundary condition cr(±oo) 2 = m 2 (2.8), we may write (A. 26) more 
compactly as 



Recall from (A. 13) that the standard wronskians W b st (k), Wf(k) vanish, respectively, 
at the bound states of H b and H c . Thus, contrary to (A.23), which holds for all 
backgrounds cr(x), we cannot have W b st (k) = Wf(k) when q — ±1, i.e., when only 
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one of the hamiltonians has an unpaired normalizable zero energy ground state. It is 
rewarding to verify that (A. 27) is consistent with these facts: 

For q = 0, H b and H c are strictly positive and thus must have the same bound 
states. Thus, we expect to find W b st (k) = W^(k), which is what (A. 27) tells us at 
q = 0. For q — 1, H b and H c are isospectral, save for the unpaired u b = ground state 
of H b . uj 2 = means k b = im. Thus, we expect to find that W b st (k) has an extra zero 
at k — im, relative to Wf(k). Indeed, from (A. 27) with q — 1 and <r(— oo) = —m, 
we find that W b st (k) = ^^Wf(k). Finally, for q = —1, the roles of H b and H c are 
reversed relative to the q — 1 case, and upon substituting q — — 1 and a{— oo) = m 
in (A.27) we find accordingly that Wf{k) = j^Wf(k). 

We note in passing that isospectrality of H b and H c is consistent with the 75 
symmetry of the system of equations in (A. 3), which relates the resolvent of D with 
that of D = — 75-D75 (recall (2.13)). Due to this symmetry, we can map the pair of 
equations (H b — u 2 )b(x, y)/cu = 5(x — y) and (H c — u 2 )c(x,y)/uj = —5(x — y) (Eqs. 
(A. 5)) on each other by 



(Note that under these reflections we also have a(x,y) <-> —d(x,y), as we can see 
from (A. 4).) The reflection a — > —a flips the signs of both asymptotic values cr(±oo) 
and thus flips the sign of the topological charge q. It changes the roles of H b and H c , 
but obviously it cannot change physics. Since this reflection interchanges b(x, y) and 
c(x, y) without affecting the physics, these two objects must have the same singular- 
ities as functions of u, consistent with isospectrality of H b and H c (save possibly an 
unpaired zero mode). 

A. 2 The Diagonal Resolvent 

Following [17, 19, 20] we define the diagonal resolvent (x \iD~ l \x) symmetrically as 



b(x,y) 



<-> —c(x, y) together with a 



—a . 



(A.28) 




(x ] — iD 1 \x ) 



C(x) D(x) J 



49 



= - lim 

2 t->o+ 



' a(x,y) + a(y,x) b(x,y) + b(y, x) 
\ c(x,y) + c(y,x) d(x,y) + d(y,x) 



(A.29) 

y=x+e 



Here A(x) through D(x) stand for the entries of the diagonal resolvent, which follow- 
ing (A. 4) and (A. 7) have the compact representation 18 

= o;6i(x)6 2 (g) = .[d x + 2a(x)} B (x) 



W b y ' 2u 

C(x) = -" Cl( ^ C2(x) , A(x)=i [d °- 2a £ )]C ( x) . (A.30) 

Clearly, from (A. 5) and (A. 6) we have 

B(x) | 1 | . 

to H b — to 2 

C(x) 



u - H^H. (A.31) 

Thus, comparing with [16], we see that B(x)/u and —C(x,uj)/uj are, respectively, the 
resolvents R_(x, to 2 ) and R + (x,u 2 ) defined in Eqs.(9) and (10) of [16]. 

The expressions for A(x) and D(x) in terms of B(x) and C(x) in (A.30) have an 
interesting property concerning zero modes. Consider, for example the case where iff, 
has a normalizable zero mode. From (A.31) we expect that as uo 2 — > 0, B(x)/uo will 
be dominated by the pole at to 2 = with a residue proportional to the wave function 
of the zero energy state (A. 18) squared, 

lim B M _ frofo) 



w 2^0+ v u 2 



(A.32) 



How does D(x) behave near that pole? Note from (A. 18) that (d x + 2a(x))b^(x) = 0. 

Thus, the operator (d x + 2a(x)), which appears in the expression of D(x) in terms 

of B(x) in (A.30), projects the zero mode b (x) out of D(x)\ Similarly, the operator 

(d x — 2<j(x)), which appears in the expression of A(x) in terms of C(x) in (A.30), 

1S A, B, C and D are obviously functions of to as well. For notational simplicity we suppress their 
explicit lo dependence. 
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projects the zero mode cq{x) out of A(x). A(x) and D(x) do not have a pole at 
uj 2 = 0, no matter what value the topological charge q takes on. 

Since B(x)/uj and — C(x, uj)/uj are, respectively, the diagonal resolvents of the 
Schrodinger operators and H c , each one of them satisfies its appropriate Gel'fand- 
Dikii [18] identity: 



-2Bd 2 x B + {d x Bf + AB 2 {(7 2 - a' - to 2 ) 



uj 2 



-2Cd 2 x C+(d x C) 2 + AC 2 (a 2 + a' - uj 2 ) = uj 2 . (A.33) 

A linearized form of these identities may be obtained by deriving once and dividing 
through by 2B and 2C. One obtains 

d 3 x B -4d x B(a 2 -a' -uj 2 )-2B(a 2 -a')' = 



dlC - A8 x C(a 2 + a' - uj 2 ) - 2C(a 2 + a')' = 0. (A.34) 
(For a simple derivation of the GD identity, see [17, 19].) 
A. 2.1 Relations Among A, B, C, and D 

We now use the supersymmetry of the Dirac operator, which we discussed in the pre- 
vious subsection, to deduce some important properties of the functions A{x) through 
D(x). 

From (A. 30) and from (A. 2) we have 

«*) = <^(-^) = ^(<^ + 

Using (A. 16) first, and then (A. 15), we rewrite this expression as 

A{x) = ^-Mi + cx&a) = ^- (hQb 2 + b 2 Qh) . 

Then, using the fact that W c = W b (Eq. (A. 23)) and (A.30), we rewrite the last 
expression as 

^_. d x + 2a f ujb 1 b 2 \ _ Xd x + 2a)B(x) 

A[x) - l ^r \w) ~ l 2^ • 
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Thus, finally, 

A(x) = D(x) . (A.35) 

Super symmetry renders the diagonal elements A and D equal. Clearly, the projection 
of zero modes (when they exist) out of A(x) and D(x), which we discussed following 
Eq. (A. 32), makes A = D possible. We certainly cannot have B and C equal, since 
when one of them has a pole at oo 2 = 0, the other does not. (See (A. 38) and (A. 39) 
below.) 

Due to (A. 30), A = D is also a first order differential equation relating B and C: 
(d x + 2a(x)) ^ = (d x - 2a(x)) ^ . (A.36) 

With the identification of the resolvents Bjuj and —C/u, respectively with R_(x, oj 2 ) 
and R + (x,u 2 ) alluded to above, we can write (A.36) as 

(2a - d x )R + = (2a + d x )R. , 

which is essentially Eq.(18) of [16]. 

We can also relate the off diagonal elements B and C to each other more directly. 
From (A. 30) and from (A. 15) we find 

After some algebra, and using (A. 23), we can rewrite this as 

-u 2 C = a 2 B + aB' + ^ 2 

The combination ujb'^/Wb appears in B" = (ujbib 2 /W b )" . After using (H b — u 2 )bi^ = 
to eliminate b" and b'^ from B" , we find 

ub^ 1 



W b 2 

Thus, finally, we have 



= L B " ~(a 2 -a' -UJ 2 )B. 



-u 2 C = -B" + aB' + (a' + uj 2 )B . (A.38) 
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In a similar manner we can prove that 

L0 2 B = --C" + aC' + (a'-oo 2 )C. (A.39) 

We can simplify (A. 38) and (A.39) further. After some algebra, and using (A. 30) we 
arrive at 

C(x) = -d x D(x)-B(x) = --^d x [(d x + 2a(x))B(x)]-B(x) 
00 Zoo 

B(x) = Ld x A{x)-C{x) = -^- 2 d x [{d x -2a{x))C{x)]-C{x) (A.40) 

(which are one and the same equation, since A(x) = D(x)). Supersymmetry, namely, 
isospectrality of H b and H C) enables us to relate the diagonal resolvents of these 
operators, B and C, to each other. 

Thus, we can use (A. 30), (A. 35) and (A.40) to eliminate three of the entries of 
the diagonal resolvent in (A. 30), in terms of the fourth. 

The relations (A.40) (or (A. 38) and (A.39)) were not discussed in [16], but one 
can verify them, for example, for the resolvents corresponding to the kink case cr(x) = 
mtanhmx (Eq. (29) in [16]), for which 

2 



C=- " , B 
1\Jm 2 — oo 2 



msechmx . 

— — - 1 



C. (A.41) 



Note that the two relations (A. 38) and (A.39) transform into each other under 

B <-> — C simultaneously with a — > —a , (A. 42) 

in consistency with (A. 28). 

From the relations in (A.40), it is clear that away from oo = 0, B and C must have 
the same singularities as functions of oo 2 . However, we see, for example from (A. 38), 
that if C has a pole l/oo 2 (which corresponds to a normalizable zero mode of H c ), 
then B will not have such a pole, and vice versa, as it should be. 
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A. 3 Bilinear Fermion Condensates and Vanishing of the Spa- 
tial Fermion Current 



Following basic principles of quantum field theory, we may write the most generic 
flavor-singlet bilinear fermion condensate in our static background as 



{^aa(t,x) T af3 ip af3 (t,x)) Icg = N J c ^ tr 



F(x\ 



CJ7 + %^ x d x — a 



x, 



reg 



( A(x) B(x) \ 
\C(x) D(x) ) 



( A B 



\ C D 



(A.43) 



where we have used (A. 29). The integration contour is specified in Section 3 (see 
Figure 1). Here a — 1, • • • , N is a flavor index, and the trace is taken over Dirac 
indices a, (3. As usual, we regularized this condensate by subtracting from it a short 
distance divergent piece embodied here by the diagonal resolvent 

(A B \ , / i<r„ 

(x\ -iD~ 1 \x) VAC = 

\ C D 



2\fm? — to 2 



VAC 



(A.44) 



V icr VAC 

of the Dirac operator in a vacuum configuration a VAC = ±m. 
In our convention for 7 matrices (1.2) we have 

/ A(x) B(x) \ 
\C(x) D(x) ) 



A{x) + D(x) Ajx) - D(x) B(x) - Cjx) B(x) + C(x) 
2 2i 7 2 7 



2 

(A.45) 

An important condensate is the expectation value of the fermion current (j M (a;)). 
In particular, consider its spatial component. In our static background <j(x), it must, 
of course, vanish identically 

(j 1 (x)) = 0. (A.46) 

Thus, substituting T = 7 1 in (A.43) and using (A.45) we find 

du 



{j\x))=iN j c £[A{x)-D{x)} 



(A.47) 
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But we have already proved that A(x) = D(x) in any static background a{x) (Eq.(A.35)). 
Thus, each frequency component of (j 1 ) vanishes separately, and (A. 46) holds identi- 
cally. It is remarkable that the supersymmetry of the Dirac operator guarantees the 
consistency of any static a(x) background. 

Expressions for other bilinear condensates may be derived in a similar manner 
(here we write the unsubtracted quantities). Thus, substituting T = 7 in (A. 43) and 
using (A. 45), (A. 35) and (A. 40), we find that the fermion density is 

, n, u » T f du r „. , _ , T r duj 2ujB(x) — id x D(x) 

(j°(x)) =iNj c — [B(x) - C{x)\ =iNj c — ^ ■ (A.48) 

Similarly, the scalar and pseudoscalar condensates are 

mx)^(x)} =Nj^ [A(x) + D(x)\ =2Nf^ D(x) , (A.49) 

and 

mxh^(x)) = N f ^ [Btx) + C{x)\ =iNd x f^RW. (A.50) 

We can also derive a simple relation between the fermion density (A.48) and the 
pseudoscalar condensate (A.50). From (A. 36) we have d x (B — C) = —2a(B + C). 
Substituting this result into (A.50) and comparing with (A.48) we deduce that 

d x (j°(x)) = -2ia(x) ($(x)>fil){x)) . (A.51) 

We may verify (A.51) easily for the kink a(x) = m tanh mx by substituting (A. 41) 
in (A.48) and (A.50). 

The result (A.51) has a simple interpretation in terms of bosonization. Since 
(A.51) stems essentially from the properties of the diagonal resolvent (A. 29) of the 
Dirac operator i0—a(x), it is independent of the 0(2N) symmetry of the GN model, 
and also independent of the dynamics of the auxiliary field a(x). Thus, we need to 
consider only bosonization of the action 

Sj, = J d 2 xip (i0 - a(x)) ip (A.52) 
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of a single Dirac fermion -0, corresponding to the Abelian case [36]. The relevant 
bosonization operator identities are 

■0-0 — > \i cos(2v / vr0) and 



i0> 7 5 0> -> /i sin(2v^F0) , (A.53) 

where is the corresponding boson, and [i is an arbitrary scale parameter. 
According to these formulas, the bosonized form of is 




(A.54) 



which leads to the equation of motion 

□0 - sin (2^0) = . (A.55) 

In the static limit, the latter equation becomes 

d 2 x <f){x) + 2^a{ x) sin (2\fK(j){x)) = , (A. 56) 

which is easily recognized as the bosonized form of (A. 51) (up to a trivial multiplica- 
tion of (A. 56) by a factor N). Thus, the identity (A. 51) is equivalent to the equation 
of motion of the bosonic field in the static limit. This of course, should be expected 
on physical grounds. After all, Eq. (A. 51) is a consequence of (A. 35) (or (A. 36)), 
which implies (j 1 (x)) = in our static background cr(x), which leads to <9o0 = 0, 
according to the second equation in (A.53). 
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Appendix B: Reflect ionless Schrodinger Operators 

and their Resolvents 



In this Appendix we gather (without derivation) some useful properties of reflection- 
less Schrodinger hamiltonians[3] and their resolvents. A particularly useful summary 
of the theory of inverse scattering in one dimension and of reflectionless Schrodinger 
operators is given in the first reference in [9] (which we partly follow here). 

Following that summary of results from the literature, we use one particular result, 
namely, the explicit formulas for a pair of fundamental independent solutions of the 
Schrodinger equation, to derive an explicit simple representation for the diagonal 
resolvent of a reflectionless Schrodinger operator with a prescribed set of bound states, 
which we have not encountered in the literature. 

B.l The Potential and Wave Functions 

Consider the Schrodinger equation 

{-d 2 x + V(x)) b(x) = k 2 b(x) (B.l) 

on the whole real axis, where V(x) is a reflectionless potential, namely, the reflection 
amplitude r(k) of (B.l) vanishes identically. Reflectionless potentials are bounded, 
and tend to an asymptotic constant value V(±oo) = Vq at an exponential rate. With 
no loss of generality we will set Vq — 0. Thus, V(x) satisfies the boundary conditions 

V(x) 0. (B.2) 

Since V(x) is bounded and tends to its asymptotic value (B.2) at an exponential 
rate, it can support only a finite number of bound states. Let us assume that the 
reflectionless V(x) has K bound states at energies 

-k\ < -«2- • • < -k 2 k < 0. (B.3) 

Eq.(B.3) is the spectrum of bound states of a one-dimensional Schrodinger operator, 
which cannot be degenerate. Thus all inequalities in (B.3) are strict. 
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The reflectionless potential V(x) is uniquely determined [3] by the asymptotic 
behavior of its K bound state wave functions only. (In order to determine a generic 
reflectionful potential, one also requires, of course, the reflection amplitude r(k).) 

A bound state wave function in a potential which tends to its asymptotic behavior 
(B.2) at an exponential rate, must also decay exponentially. Thus, from (B.l), the 
n-th bound state wave function if) n (x), which corresponds to bound state energy — k^, 
decays asymptotically as 

i> n {x) — > = c n exp-K n x, (B.4) 

with some real coefficient c n . The parameter c n is determined from the requirement 
that ip n be normalized to unity. With no loss of generality we take c n to be positive. 
For convenience, let us name the right hand side of (B.4) as 

X n (x) = c n exp -n n x . (B.5) 

The asymptotic behavior of ip n (x) is determined by two positive parameters K n 
and c„, and thus, there are 2K parameters at our disposal: 

Ki > k 2 > • ■ ■ kk > and ci, c 2 , • • • ck ■ (B-6) 

These are the "scattering data" alluded to in Section 3 (see (3.4)). These parameters 
determine the reflectionless potential V(x) with K bound states (B.3) neatly as 

V(x) = -2— logdetA, (B.7) 
where A(x) is the symmetric K x K matrix 

A mn = S mn H — . (B-8) 
It can be shown that the K dimensional vectors ip n and \ n are related by 

K 

^A mn i) n = \ rn . (B.9) 

n=l 

We can solve this linear equation and thus obtain an explicit formula for the if) n (x). 
The desired solution is 

, . . 1 detA^ , n . 
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where A^ n \x) is the K x K matrix obtained from (B.8), by replacing the n-th column 
of A by its derivative. 

Inverse scattering theory also provides us with explicit expressions for the pair 
of fundamental independent solutions of (B.l). Let {bi(x) , b 2 (x)} be a pair of in- 
dependent fundamental solutions of (B.l) subjected to the (standard, or canonical) 
asymptotic boundary conditions 



-ikx 



b 2 (x) — > • 



ikx 



;b.ii) 



(i.e., Eq. (A.9) with A\ 



(i) 



I (2) 



1). Then, inverse scattering theory tells us that 



for ImA; > 0, b\ and b 2 are given by 19 

1 



biix) 



ikx 



t(k) 

b 2 (x) = e lkx 



1 , . ^ \ n (x)i) n (x) 

n- 



where 



1 _ . J2 ^n(x)jj n (x) 



b 2 {x)b' 1 (x) - bi(x)b' 2 (x) 



(B.12) 



(B.13) 



t(k) 2k 

is the inverse of the transmission amplitude t(k) of (B.l). It is proportional to the 
wronskian 

W = b 2 {x)b\{x) - b^b'^x) (B.14) 
of b\ and b 2 , as was mentioned in Appendix A, 



1 



iW 



(B.15) 



t(k) 2k ' 

According to inverse scattering theory, the transmission amplitude t(k) of the 
reflectionless potential with K bound states which correspond to the data (B.6) is 
simply 

K h 4- itc 

(B.16) 



t (k) = n 1 f ;/ " 



_2 k 



19 The functions b\ and 62 are, respectively, the functions 4>2 and (f)\ in the first reference cited in 



[9]. 
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This transmission amplitude is a pure phase, as it should be, since the potential V(x) 
is reflect ionless. Note the simple poles of t(k) at the bound states, or equivalently, 
the simple zeros of l/t(k), which arise, as we have already discussed (in connection 
with (A. 13)), since l/t(k) is proportional to the wronskian of the two fundamental 
solutions, and both these solutions become proportional to the bound state wave 
function, and thus linear dependent. 

Thus, as k tends to one of its bound state values, say at k — iK m , the sum in the 
square brackets in the expression for bi(x) in (B.12) has a pole, which dominates the 
sum . But at the same time, l/t(k) vanishes linearly as k — in m (up to some constant 
one can calculate from (B.16)). This is, of course, a general feature of scattering 
theory, and is not particular to reflect ionless potentials. Thus, at the limit k — > iK m 
the expression for bi(x) in (B.12) collapses simply to 

bi(x) — const. ip m (x) (B-17) 

(where we used e KmX X m (x) = c m ). To see that at k — in m , the other fundamental 
solution b2(x) in (B.12) is also proportional to ip m (x) already requires the special 
properties of reflectionless potentials. Indeed, at k — iK m , the expression for b2(x) 
may be written 



bo (x) = e 



Ft. f-,1 X 



1 K 

T E (A mn — 5 mn ) 1p n (x) 
Am n=1 



— i) m {x) , (B.18) 

Cm 



where we have used (B.5) and (B.9). bi(x) and 6 2 (x) are indeed both proportional to 
the bound state wave function at k — iK m . 

We can also use the explicit forms of the fundamental solutions (B.12) and of the 
transmission amplitude (B.16) to derive an explicit formula for the integral over the 
potential V(x). According to (B.7), the definition of the matrix A^ n \ and (B.10), we 
have 

J V{x) dx = 2 **% } = "2 lim E Ux)M*) ■ (B.19) 



71=1 



Here we used the fact that lim A n (x)ib n (x) = 0, obviously. However, lim \ n (x)ib n (x) 
is a constant, which has to be determined. In fact, in (B.19) we need only the sum 
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of these constants. This sum can be extracted from (B.12). Indeed, from (B.ll) and 
from (B.12), we observe that 



/ K 

t(k) = 1 + % hm PT 

rr^-oo \ n =l 



Thus, combining (B.20) and (B.19), we finally arrive at the sum-rule 

dk 



(B.20) 



/ V(x) dx = 2i I ^-.t(k) , 
J Jc 2m 



fB.21) 



where the contour C encircles all the poles of t(k) in (B.16) in the upper half /c-plane. 



The Diagonal Resolvent 
We now use (B.12) to derive a simple formula (Eq. (B.28) below) for the Green's 
function of the reflectionless hamiltonian 



H = -d 2 x + V(x 



(B.22) 



in (B.l). We have not encountered the formula we are about to derive in the literature. 

In a similar manner used to derive the Green's functions in (A. 7), we may express 
the Green's function of H in terms of the fundamental solutions (B.12) as 



r(<r „■ h\ - M 1 u,\ - 0(x-y)b 2 (x)b 1 (y) +9(y-x) b 2 {y)h{x) 
u{x,y,K) — \x\— —\y) — 



H-k 2 

Then the diagonal resolvent of H, 

R(x;k 2 ) = G{x,x;k) = 

(compare with Eq. (A. 29)) is 
bi(x)b 2 (x) i 



W 



(B.23) 



lim G(x, x + e; k) + G(x + e, x; k) 



R(x; k) 



W 



2k 



K \ nh 
1+^ 



7 : — 

m=l k + %Kr - 



(B.24) 



~t\ k iKn_ 

The last expression may be simplified considerably by simple algebraic manipula- 
tions. First, let us multiply the factors in (B.24). Thus, 



R(x; k) 



2k 



K 



n= i \k iftn k + 



1 \ K 

) A n ^n + 

' lK n/ m,n=l 



(k - tK n )(k + 



(B.25) 
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Then, concentrate on the last term in (B.25). Using (B.8), we can write it as 

(k-in Mk + iK ) = ~ l{Amn ~ 6mn) If^T " ITuT ) ^ n ■ (R26) 

Then, summing the last expression over m,n and using (B.9) we obtain 

^m^n^Pm^Pn _ ^ ^ / 1 ^ A ^ 2 ^ Kra ^ 

m n= i — if^n)(k + ift m ) n=1 \/c — i/t n A; + in n J n=1 A; 2 + /t 2 

(B.27) 

Finally, substituting (B.27) in (B.25) we arrive at the desired simpler formula for the 
diagonal resolvent: 

^^(l-SSjgL). (B.28) 

As one trivial consistency check on (B.28), note from (B.28) that R(x; k), as a function 
of energy E = k 2 , has a simple pole at each of the bound state energies E n = —n 2 n) 
with residue — ip 2 , as it should be. 20 

Also note that (B.28), being the diagonal resolvent of (B.22) with potential (B.7), 
must satisfy the GD identity 

-2RR" + R' 2 + AR 2 {V{x) - k 2 ) = 1 . (B.29) 

In fact, in [16], the analysis went the other way around, and used (B.29) to derive 
(B.28) in the case of a single bound state, K — 1. 

As yet another consistency check on (B.28), note that as \x\ — > oo, it tends to 
R = i/2k, which is indeed the asymptotic solution of (B.29) in view of (B.2). 

B.2 Application to the Dirac Operator 

In order to make contact with the discussion in the text, let us identify the hamiltonian 

(B.22) with one of the operators H b or H c in (A. 6). With no loss of generality we 

chose to identify it with H b = —d 2 + a(x) 2 — a'(x). Actually, due to the boundary 

20 Note that (B.28) is not the most obvious way to write a function of x and k 2 with these properties 
(e.g., the function similar to (B.28), but with all prefactors iHn/Vh^ removed, has these properties, 
and also the correct large k behavior i/2k, and is much simpler). Thus, we really had to go through 
all the steps of the derivation to obtain (B.28). 
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conditions (2.8), the potential a(x) 2 — a'(x) tends asymptotically to m 2 and not to 
zero. Thus, in order to comply with (B.2), we consider H = H b — m 2 , so that the 
eigenvalue equation (A. 8) H b b = uj 2 b turns into (H b — m 2 )b = (uo 2 — m 2 )b, which 
coincides with (B.l), due to oo 2 = k 2 + m 2 (A. 10). At the bound state with k 2 = —k^ 
we thus have u 2 = m 2 — k^. Thus, from (B.3) the bound state energies of H b are 

< J* < toj ■ ■ ■ < to 2 K < m 2 . (B.30) 

Since (B.3) is non-degenerate, all the u 2 must be different from each other. 

In the topologically trivial sector of the model, q — 0, H b — —d 2 + a(x) 2 — 
a'(x) and H c = —d 2 + a(x) 2 + a'(x) are strictly isospectral, and also reflectionless. 
Thus, we have two different reflectionless potentials, with the same set of bound state 
energies (B.30). Clearly, it is the additional K parameters ci,---ck that remove 
the ambiguity between these two isospectral potentials. The parameters c±, ■ ■ ■ ,ck 
parametrize a K dimensional family of reflectionless potentials with given K bound 
state energies (B.30). As the c n vary, the system wanders around in that family 
space without changing its bound state energies. In particular, H b and H c should 
correspond to two points in the K dimensional space of the c n 's. Let us denote 
these points, respectively, as {c n (H b )} and {c n (H c )}. It will be interesting to find a 
representation of the supersymmetric transformation H b <-> H c as an invertible map 
on the K dimensional vectors 

{c n (H b )} <- {c n (H c )} . (B.31) 

Computation of this transformation in the simplest case, i.e., when there is a single 
bound state, K = 1, is given in section B.3 (see (B.47)). 
Finally, from (A. 31) we have 

^-Wa^W-lHS^)- (B ' 32) 

We can then use (A. 38) to determine C(x). (Alternatively, due to (A. 31) and isospec- 
trality of H b and H c , we can just write an expression for —C(x) similar to (B.32) with 
ip n being the bound state wave functions of H c .) 
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B.2.1 The Condensate a(x) 



Let us now determine the condensate a(x). We will consider the cases of zero topo- 
logical charge and non-zero topological charge separately. 

(a) topologically trivial condensates, q = 
In this case, none of the operators H b and H c has a normalizable zero mode. Thus, 
none of the u; n 's in (B.30) vanishes. ( Equivalent ly, none of the re n 's in (B.3) can be 
equal to m.) Since Hb does not have a normalizable zero mode, the fundamental 
solutions (B.12) at oo 2 = 0, i.e., at k — im, 



b i( x )\ k=im = 



b 2 {x)\ 



1 



t(im) 



1 + V ^ n ( x )^ n ( x ) 



n=l 



m — K r 



^ \ n {x)i) n {x) 



(B.33) 



are not square-integrable functions. They are the two independent solutions of 
Q^Qh,2 = 0. According to (A. 18), one of these solutions must be b (x) = exp — j x a{y)dy. 
Thus, if we knew bo, we could determine 



cr(x) = --^\ogb (x) . 



Clearly, 

On the other hand, 
bi(x)\ k 



b (x) — > e 



and 



-cr(oo)x 



(B.34) 
(B.35) 

(B.36) 



'1-.+00 t{im) 

since in that limit, the terms proportional to \ n (x)ip n (x) in (B.33) are negligible 
compared to 1. Now, identification of bo(x) with one of the functions in (B.33) is just a 
matter of comparing (B.35) and (B.36). Thus, if a (00) = —m, then b (x) = bi(x)\ k=im , 
and if <r(oo) = m, then b (x) = 6 2 (^)| fc=im - Therefore, we conclude that 



a(x) = —m — — log 
ax 



a(x) = m — — log 
ax 



1 + J2 ^n{x)i>n{x) 



n=l 
K 



m — K r 



\ n {x)ip n (x) 



a 00 = — m 



^ m + K n 



cr(oo) = m , (B.37) 
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or more compactly, 

a(x) = a(oo) — — log 
dx 

This is manifestly a configuration belonging to the q = sector, since the expression 
on the right hand side of (B.38) tends to er(oo) as x — > ±00, i.e., on both sides of the 
one dimensional world. 



^ _ \ n (x)ip n (x) 

hi K n + CT(OO) 



(B.38) 



(b) the q = 1 sector 

According to the discussion following Eqs. (A. 18), (A. 19), in the q — 1 sector, only Hj, 
has a normalizable zero mode. According to (A. 18) this normalizable ground state is 
bo(x) = exp — J x a(y)dy, and thus must be just the lowest energy state tpi in (B.10), 
corresponding to uj\ = 0, i.e., K\ = m. Thus, from (B.10) and (B.34) we find 

.(,) = --, og* lW = -- log (--— j . (B.39) 

Let us make the side remark, that due to the discussion which led to (B.17) and 
(B.18), the present case can be considered as a singular limit of the previous case (a) 
(with er(oo) = m, as appropriate to the sector q — 1), simply by taking the second 
equation in (B.37) at K\ — > m. Following the same steps as in the derivation of (B.18) 
with n m = K\ = m we end up with 

a(x) = m - Ai og ^ = _i_ tog^(x) (B.40) 
dx Ai dx 

i.e., back to (B.39). 

(c) the q — — 1 sector 

According to the discussion following Eqs. (A. 18), (A. 19), in the q — — 1 sector, only 
i7 c has a normalizable zero mode. If it has K — 1 additional bound states at positive 
energies, they will be isospectral with those of H b . Thus, in that case H c has K bound 
states at energies 

= u\ < (J*--- < u 2 K < m 2 . (B.41) 
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Determination of a a(x) configuration that carries topological charge q — — 1 can be 
reduced to the previous case, of determination of the configuration — a(x) that carries 
q — 1. A sign flip cr(x) — > —cr(x) obviously interchanges Hf, and H c . Thus, —cr(x) will 
give rise to an with bound state spectrum (B.41), and we will be able to determine 
this —cr(x) by the method of case (b), from (B.39). 

So far, we have assumed K > 2 in (B.41). If K — 1, the zero mode is the only 
bound state of H c , and cannot have bound states at all. In that case, it is therefore 
the free particle hamiltonian H b = —d 2 + m 2 . Thus, a 2 — a' = m 2 and a(x) is the 
CCGZ antikink a(x) = — mtanhmx. 

B.3 Reflect ionless Potentials with One and Two Bound States 

We end this appendix by working out the cases of reflectionless potentials with one 
and two bound states, which are the two cases relevant to our discussion of stable 
static bags in the GN model. 

B.3.1 A Single Bound State, K = 1 

Let us set the bound state energy at uj 2 < m 2 (with the corresponding k 2 = m 2 —lo 2 ). 
We also have a single parameter c > 0. Then, from (B.5), (B.8), (B.10) and (B.7) we 
find 

2 — Ikx 

A(x) = ce- ra , A(x) = l + ^- (B.42) 

and 

MX) = ~l = /I S6Ch {<X ~ Xo)) (R43) 

where xo is determined from 

e~ KX0 = \[2kjc. (B.44) 

■ipb is of course normalized to 1. From (B.7) we find the potential a 2 — a' = m 2 +V(x) = 
m 2 — 2d 2 log A(x) = m 2 — 2k 2 sech 2 (k(x — x )), which is a potential of the Poschl- 
Teller type. 

The role of the parameter c is explicit in these formulas: it just shifts the center 
of the potential, and thus cannot affect the energy of the bound state. Thus, log c is 
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essentially the translational collective mode of the static soliton a(x), which is obvious 
from the explicit expression (B.46) below. 

The diagonal resolvent is found from (B.32) and is given by 

2 sech 2 (k(x — Xq)) 



K J 2k I u 2 - Ji 



ioj 
2k 



K 



uj 2 -Ji 



(B.45) 



in agreement with the expression found in [16]. 

Let us determine the corresponding a(x). Since u b > 0, we have to use (B.38). 
We find 



a(x) = <j(oo) — log 

ax 



1 - 



X(x)ip b (x) 



K + (j(oo) 

cr(oo) + /ttanh [k(x — xq)] — /ttanh 



k(x-X ) + - log 



/ m + 


")] 


Vm — 


kJ . 



(B.46) 



in agreement with the result quoted in Eq. (3.28) of [1]. It has the profile of a bound 
state of a kink and an antikink, with interkink distance ^- log (s2±s.) m It thus carries 

' 2k ° \m—K J 

topological charge q = 0. 

Let us now find the representation of the supersymmetric transformation <-> 
H c , alluded to in the previous subsection, as a map (B.31) between the two nor- 
malization constants c(Hb) and c(H c ) in (B.31). Since the transformation <-> H c 
is achieved simply by cr(x) <-> — cr(x), we have to find a pair of normalization con- 
stants c(H b ) and c(H c ), which will yield, upon substitution into (B.38) and (B.42) 
(with boundary conditions er(oo; b) = m = — cr(oo; c) in the latter equation), two a(x) 
configurations of opposite signs. The calculation is straightforward, and we find 



c\H c ) = c\H b ) 



m + k 



or equivalent ly, from (B.44), 



kx (H c ) = Kx (H b ) + ]- log m 

I m + k 



(B.47) 



(B.48) 



It is easy to verify, by substituting (B.48) in (B.46), followed by cr(oo) = m — > —m, 
that a(x — xo(Hb)) = —a(x — xq(H c )), as required. Furthermore, note that if xq = 
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x (H b ) in the first hyperbolic tangent in (B.46), the second hyperbolic tangent there 
is shifted by x (H c ). 

As k tends to m, u 2 tends to zero. An eigenvalue u; 2 = cannot occur in the 
topologically trivial sector. Thus, the limit K->min (B.46) is singular, and the kink 
and antikink become infinitely separated, as we discussed following (3.36). 

In the sector with topological charge q — 1, we have u 2 = in the spectrum of 
H b , and we should use (B.39) for a(x). Thus, from (B.39) and (B.43) with k = m we 
find 

cr(x) = — — log^i(x) = mtanh (m(x — xq)) , (B.49) 

which is of course the CCGZ kink. In the CCGZ kink background, Hb has a single 
bound state at Ub = and no other bound states. Thus, H c has no bound states at 
all, and is thus the free particle hamiltonian H c = —d 2 +m 2 . Indeed, with (B.49), one 
has a 2 +a' = m 2 for the potential of H c . The CCGZ antikink was already mentioned 
at the end of case (c) in the previous subsection. 



B.3.2 Two Bound States, K = 2 



We now determine the reflectionless potentials with two bound states. Let us set 
the bound state energies (B.30) at < uo\ < uj\ < m 2 , with the corresponding 
< k 2 < Ki < m. We also need ci, c 2 , the two positive parameters in (B.4). The case 
relevant for our discussion of stable static fermion bags with two bound states is, of 
course, the case = uo\ < u>2- We will discuss it at the end of this subsection. 
Then, from (B.5) and (B.8) we construct 



We will also need 



c\ e 



-2k\x 



( 1 _|_ c i c -2k\x 
L T 2ki 



C1C2 c -(ki+K9.)x 



_ClC2_ e -(K 1 +K 2 )x \ 
K1 + K2 



_Cl£2_ ( 

K1+K2 



-{k\+K2)x \ 



1 + ^e~ 2K2X 



(B.50) 



y -cic 2 e - ^ 1 ^ 2 ^ l + ^e- 2K2:r 



1 + ^ 2K1X 



_£l£2_ g— (ki +K2 )x 

K1+K2 



-(ki+k,2)x \ 



-cic 2 e 



„2 „—2k2X 

c 2 e 



;b.5i) 
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The determinants of these matrices are 

2 2 2 2 2 

F(x) = det A = 1 + ^ e - 2 ^ x + ^ e ~ 2K2X + ( ^U^.) e -2(«i+« 2 )x 

2«i 2k 2 4kik 2 Vki + K2J 

F 1 (x) = detA^ = -cle- 2 ^(l + ^^^e- 2K A 

F 2 (x) = Aet = -c 2 2 e~ 2 ^ x (l + ^^^e- 2K A . (B.52) 
V ; 2 V 2kik 2 + ki ; V ; 

The potential cr 2 — cr' can be found from (B.7) and (B.52) as 

d 2 

a 2 - a' = m 2 + V(x) = m 2 - 2— logF(x) . (B.53) 

ax 2 

We will not need an explicit expression for V(x). 

The bound state wave functions may be found from (B.10) 

1 F 1 (x) 



X 1 (x) F(x) 

( 1 + ° 2 Kl 

\ 2ft 2 Kl 



3 K1X I 1 _(_ c 2 Kl— K2 „ — 2k<2X 
1 > ^ ... -i-K2 



I i c i c -2k,i.t _i_ _52_g-2K2a; _|_ c i L 2 ( ki-k 2 j g-2(Ki+K2)a: 
2ki 2k 2 4kiK2 \ki+K2/ 

1 F 2 (aQ 

A 2 (x) 

>e~ K2X (l + ^Lt^i e ^ lX \ 

(B.54) 



]_ _|_ _£i_g-2/cia; _i_ _£Lg-2K2a; _|_ c i c 2 / ki-k 2 ~\ 2 g-2(Ki+K 2 )z 
2ki 2k 2 4kiK2 \ki+k 2 / 

The parameters ci, c 2 appear in the potential and wave functions (B.53) and (B.54) 
explicitly. They parametrize a 2 dimensional family of reflect ionless potentials with 
given two bound state energies uj\ < uj\. In other words, as C\ and c 2 vary, the 
profiles of V(x), ip±(x) and ^2 0*0 (and in particular, their "center of gravity") change 
accordingly, with the energies fixed. Thus, as in the single bound state case, C\ and c 2 
represent collective coordinates of the soliton a(x). This is manifest, for example, in 
the explicit expression (B.64) below, for the a(x) corresponding to K\ = m > k 2 = k. 
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One particular choice of c\ and c 2 leads to a symmetric potential V(— x) = V(x). 
From (B.7) and (B.52) we see that this will occur when log F f[x) = ^' *- e -> a ^ 
F(—x) = e ax+b F(x) for some constants a, b. Setting x = in the last equation tells 
us that 6 = 0. Thus, an even potential occurs when 

F(-x) = e ax F(x) . (B.55) 

Making this demand on F(x) in (B.52), leads, after some algebra, to the condition 

_4_ = _cf_ = (B 56) 

For these values of C\ and c 2 we obtain from (B.52) 

F(x) = 1 + Kl + K2 ( e -^ix + e ~ 2K2X ) + e -2(K 1+ K 2 )x _ ( R57 ) 

It is straightforward to check that (B.57) satisfies (B.55) with a = 2(k± + k 2 ), 

F(-x) = e 2(Kl+K2)x F(x) . (B.58) 
In this case, we can read the ground state wave function off (B.54). It is 

fr(*) = - ^ (B-59) 

Similarly, the excited state is 

*^ £2 ^4m — -■ (R60) 

Thanks to (B.58), ip\{— x) = ipi(x) is symmetric, and ip2(— x ) = —^(x) is antisym- 
metric, as it should be. 

Avoiding Degeneracy It is amusing to observe how this one-dimensional quantum 
system avoids degeneracy in its bound state spectrum, as it must. (From now on, we 
consider generic backgrounds and do not assume a symmetric V(x).) If we attempt 
to construct a reflectionless potential with two degenerate bound states, i.e., with 
Ki = k 2 = k, we observe from (B.54) that 

-01 (x) 4>2{x) e~ KX 



Ci Co 1 i c j+ c 2 c -2kx 
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(B.61) 



The two wave functions become linearly dependent, and thus, there is no degeneracy. 
Rather, the system loses one bound state, since now 

2 2 

F(x) = 1 + c -l±^l e -^ (b.62) 

coincides with A(x) = 1 + c e 2 ™ in (B.42) of the single bound state case, with 
c 2 = c\ + c\. 



The Case = uj\ < uj\ Rather than continuing the general discussion, let 
us specialize at this point in the reflectionless potential which is relevant for our 
discussion of stable fermion bags with two bound states, namely, the case with = 
uj\ < lo\. This corresponds to setting K\ = m > k 2 . Let us also denote k 2 = k in 
(B.50) through (B.54) above. In particular, the ground state wave function is 



A (x) = 5 ^ J - 2 . (B.63) 

-y i £^ e -2mx _i_ 22.p-2.KX _i_ c i c 2 / m-K \ „-2{m+n)x 
2m 2k ArriK \m+K J 

The wave function tpi(x) is the normalizable zero mode of H b . Thus, we are in the 
q — 1 topological sector. In this sector we can determine a(x) from (B.39). Thus, 
substituting (B.63) in (B.39), we obtain 

a(x) = ——\ogipi(x)=m + 



dx 1 l+ 2nrn+K e 2KX 



cS m—K 



2m,K m+K g2mx _i_ 2mK rn+K ^2kx 



_ / \ C. [Tfl — K)~ < \ in — 'i 

- 2(m + K) g 2 — 2 • 

I I 2m / jtj+kA 2ml i 2k / m+ K \ 2kx i 4r»K / m+« \ 2(m+K)a; 
c| \m—K J c| \m—K J c'fci, \m—Kj 

(B.64) 

One can readily check that the asymptotic values of a(x) are cr(oo) = m and cr(— oo) = 
m + 2k — 2(/c + to) = —to. Thus, a(x) indeed carries topological charge q — 1. 

By varying the parameters Ci,c 2 (and keeping k fixed), we can modify its shape 
(e.g., translate its center of gravity), without affecting the bound state energies. Eq. 
(B.64) represents a two parameter family of isospectral kink backgrounds. Thus, as 
we have asserted earlier, c\ and c 2 are related to the collective coordinates of the 
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soliton a(x). In this case, we have two such translational collective coordinates, xo 
and y , given by 

e 2mx = and e 2 Kyo = £2. (Rg5) 

2m 2k' v ; 

in complete analogy with (B.44). In terms of x Q and y , we can write (B.64) in a 
slightly less cluttered form as 

2k 

<7{X) - m + m±K e 2 K ( x - yo ) 
m—K 

— 2 m + AC — — — — - 



m+K, \ e 2m(x-x ) i ( m+K. \ e 2 K (x-y Q ) i ( m+« \ g2m(x-x )+2«(x-j/o) 

(B.66) 

We may simplify (B.63) and (B.66) further as follows. In terms of the distance 
scale 21 



2k \m — k 
we rewrite (B.63) as 

/— cosh - y + -R)] 

u = V2m 



e Ki? cosh [m(x — x ) + k(x — y ) + 2kR] + e nR cosh [m(x — x ) — k(x — y )] 

(B.68) 

Thus, we obtain the somewhat more transparent expression 
a(x) = — — \ogipi(x) = — k tanh [k(x — yo + R)] 

(JjJL 

sinh [m(x — x ) + k(x — yo) + 2kR] + sinh [m(x — xo) — k(x — yo)] 
b e -nR pQg^ [ m ( x _ _|_ K (x _ _|_ 2kR] + e KR cosh [m(x — x ) — k(x — yo)] 

(B.69) 

for cr(x), where we used (m + K)e~ KR = (m — /t)e Kfi = \Jm 2 — k? = cut,. 

With no loss of generality, we can set, of course, one of the collective coordinates 
to zero. In [23] we have studied the profile of (B.66) in some detail. (See Figures 1-3 
in [23].) 

21 Note that (B.67) coincides with the interkink distance in (B.46). 
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We can also tune the parameters c\ and C2 in (B.64) (for purely mathematical 
reasons and without apparent physical motivation), such that Vb(x) = a 2 (x) — a\x) 
be an even function (and a(x) be an odd function). Thus, using (B.56) and (B.65), 
we set 2mx = 2ki/ = 2kR in (B.69), and obtain 

sinh Urn + k)x] + sinh Urn — k)x] 
alx) = — K,tanh{KX) + cu b 5 — r r-^ B —r-. --r . (B.70) 

e -KH gQg^ _|_ K j^j _|_ 6 kh gQg^ _ K j x ^ 

Then, if for example, we set k = m/2 in (B.70) we see that it simplifies considerably, 
and we obtain 

a (x) = mtanh (mx/2) , (B-71) 

i.e., a deformation of the CCGZ kink mtanh mx. It is clearly spread out in space as 
twice as much as the CCGZ kink. 



73 



Appendix C: Reflectionless cr(x)'s Are Extremal 



As a consistency check of our calculations in Section 3, we will verify now that van- 
ishing of the right hand side of (3.19) guarantees that the corresponding a(x) is 
indeed extremal among all possible admissible static configurations (i.e., a solution of 
J|j|y = 0), and not just an extremum over the space of reflectionless configurations: 
Substituting V b and V c from (3.9) into (3.10), and integrating by parts over 5a'(x), 
we obtain 



58[a} = Jdx5a(x){^ + N 



g 2 J 2vr 



i(2a(x) + d x )B(x) i(-2a(x) + d x )C(x) 



2u 2u 

(C.l) 

With the help of (A. 30) we recognize the two expressions in the square brackets 
in (C.l), respectively, as D(x) and A(x), the two diagonal entries of the diagonal 
resolvent (2.11). Thus, from (C.l) we conclude that 

58 [a] a(x) f duj 



+ N J £ (A(x) + £>(*)) . (C.2) 



5a(x) g 2 

Finally, we may simplify (C.2) further, by invoking (A. 35), which tells us that A(x) = 
D(x). Thus, we may write the variational derivative as 

58\ a] a(x) f du 



2N I £ D{X) ■ (C - 3) 



5a(x) g 2 

As a consistency check, we obtain (C.3) in the next subsection as the static limit of 
the general extremum condition (2.3). 

For later use, let us record here (C.3), evaluated at the vacuum condensate a = 
(J VAC = ±m. Thus, 

^+2^/^ = 0, (C.4) 

where D VAC = 2 ^ 2 AC ^ 2 is given in (A. 44) in Appendix A. It is a trivial matter to 
verify from (A. 44) and (2.3) that (C.4) is equivalent to the bare gap equation (2.4). 

Equations (C.l) - (C.3) are general identities, valid for any static a(x). We will 
now evaluate (C.3) at a reflectionless a(x) and show that its right hand side vanishes. 
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We start by substituting the relation D(x) = i &±M?MM (f rom (A.30)), and the 
expression (B.32) for B(x) into (C.3), and obtain 



The first term on the right hand side of this equation vanishes due to the gap equation 
(C.4) (compare with (3.16)). As for the sum in (C.5), we now show that it vanishes 
term by term. One has to be a little bit more careful in case H b has its ground state 
at zero energy. Consider first all terms in the sum which correspond to positive bound 
state energies cj 2 > 0. Due to (3.19), each one of these terms is proportional to the 
corresponding J^, which vanishes at an extremal uj n . If has a zero energy ground 
state, i.e., if oj\ = 0, its contribution to the sum in (C.5) vanishes also, but due to 
the fact that it is proportional to (d x + 2a(x))ipf. The latter expression vanishes, 
since in this case ipi(x) is proportional to b (x) = exp — J x a{y)dy in (A. 18). In any 
case, the right hand side of (C.5) vanishes. Thus, we have verified that a reflectionless 
<j(x; u n , c n ), with parameters u n that satisfy (3.19) (and any set of c n 's), is an extremal 
configuration of £[a\. 

C.l A Consistency Check of the Static Saddle Point Equa- 



As a consistency check on (C.3), let us verify that the saddle point equation -^j- = 0, 

with the variational derivative given by (C.3), is equivalent to the extremum condition 

(2.3), evaluated at a static configuration 22 . 

To this end, observe that for static a(x) backgrounds, we have the (divergent) 

22 Note that (2.3) is obtained by taking a generic space-time dependent variation 5<r(x, t) of the 
effective action S e fj around a static <r(x), whereas (C.3) is obtained from varying the energy func- 
tional, which is defined, of course, only over the space of static <r(x) configurations, and thus allows 
only static variations 5a(x). 





(C.5) 



tion (C.3) 
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formal relation 



1 r duj . . 

M^-M = J-( X \ 



1 



\x) 



. f duj 



(C.6) 



CJ7 + i^ l d x — a 
( A(x) B(x) \ 
\C(x) D(x) ) 

(see Eq. (A. 43)). Thus, we may write the (bare) saddle point equation (2.3) for static 
bags as 



a{x 



N J^(Mx)+D(x)) = (C.7) 



5a(x,t), q 2 

v ' ' | static ° 

which is manifestly equivalent to the saddle point condition obtained by equating the 
right hand side of (C.2) to zero. 

From (A. 43) (with T = 1), we observe that 



N 



duj 
2^ 



(A(x) + D(x)) = N 



2tt 



tr 



C(x) D(x) 



(tp a (t,x)ip a (t,x)) la(x) (C.t 



Thus, (C.7) is just the statement that 

{Mt,x)ip a (t,x)) la(x) 



a (x) 

IF 



(C.9) 



i.e., that the auxiliary field a in (1.1) is proportional to the condensate (ipip), which 
is nothing but the equation of motion of a with respect to the action S in (1.1). 
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Appendix D: The 0(2N) Energy Multiplets 



D.l Quantization of Majorana Fields in a Background a(x) 

The 0(2N) quantum numbers of the bound state multiplets associated with the static 
solitons discussed in Sections 3 and 4 can be determined by considering the action 23 

S = J d 2 x\^ti, a (i$-a)^ (D.l) 

for N non-interacting Dirac fermions in a static background a(x). Following (1.4), 
we write (D.l) in the explicit 0(2N) invariant form 

- 2N 

S= d 2 xY: [i&ti - i&<ri& ~ (gnZMx)] , (D-2) 
J i=i 

in terms of the 2N Majorana fermions where we used the obvious notations 

<Pa = 6a-l, Xa = 6a (D.3) 

in (1.3). 

Quanta created by the Dirac-Majorana fields £j may form non-interacting, many- 
fermion bound states in the external background cr(x), which fall into multiplets of 
the 0(2N) symmetry group. In this appendix we study these multiplets. 

In the context of the large-iV limit of the GN model, a(x) is determined self- 
consistently from the static saddle point equation (C.7). The 0(2N) multiplet content 
associated with that soliton then follows from the simple group theoretic considera- 
tions presented here. 

The Dirac-Majorana equation associated with (D.2) is 

i£-i<Ti?-(<T20<r(x)=0. (D.4) 

Stationary solutions f n (x,t) = e~ %u)nt f n (x) of (D.4) with uj n ^ 0, are complex. 

However, since (D.4) has purely imaginary coefficients, the spinors Ref n (x, t) and 

23 The presentation in subsections D.l and D.2.1 follows, in part, portions of section 3 of [30]. For 
more details on the representation theory of orthogonal groups see, e.g., [37]. 
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haf n (x,t) comprise a pair of independent solutions of (D.4). This is just the state- 
ment, in terms of Majorana spinors, that the stationary solutions of the Dirac equation 
(1.6), [i0 — a(x)]ip = 0, come in charge-conjugate pairs e~ luJnt f n (x) and e luJnt f*{x), 
as was discussed in the Introduction. In addition, if a(x) is topologically non-trivial, 
as we discussed in subsection A. 1.1, (D.4) has a non-degenerate real zero energy 
eigenstate fo(x) [31]. Thus, we need only consider stationary solutions of (D.4) with 
u) > 0. 

From (D.2) we see that the quantum Dirac-Majorana field operators t),i = 
1, . . . 2N satisfy the equal time canonical anticommutation relations 

{&a(x, t),Zjp(x', t)}= ^SijS a pS(x - x') , (D.5) 

where a,f3 are spinor indices. As usual, these (hermitean) operators are to be ex- 
panded [30] in normal modes of (D.4). This expansion, which is exact for the simple 
field theory (D.l), is the semiclassical expansion of the fermion field in a self-consistent 
static background in the GN model. For example, in the presence of a topologically 
non-trivial cr(x), with its unpaired real zero energy bound state spinor fo(x), we have 
the expansion 24 

£ia(x,t) = foa(x)bi + {fna(x,t)a nji + /* a 0M)4,;} • (D.6) 

w n >0 

(For topologically trivial backgrounds, we will obtain an expansion similar to (D.6), 
with the first term fo a ( x )bi excluded.) 

In (D.6) a n i and a) n j are, respectively, annihilation and creation operators for a 
particle of type % in the state corresponding to f n (x,t). Indeed, from (D.5) (with the 
appropriate normalization of the complete set of stationary states f n ) one obtains the 
anticommutation relations 

and all other anticommutators involving a's or a^'s vanish. 

24 For simplicity, we assume here that the system lives in a big spatial box, rendering its spectrum 
discrete. 
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The operators 6, are hermitean and transform according to the vector representa- 
tion of 0(2N), as they inherit these properties from the Moreover, the canonical 
anticommutation relations (D.5) (with the appropriate normalization of /o) imply 
that the 6j satisfy the Clifford algebra 

{b i ,b j } = 26 i , j . (D.8) 

In other words, the bi have the transformation properties and the anticommutation 
relations of the gamma matrices of 0(2N). 

In terms of the creation and annihilation operators discussed above, the (normal 
ordered) Hamiltonian corresponding to (D.2) is 

H = J2 w » fl U,i' ( D - 9 ) 

D.2 The 0(2N) Multiplets 
D.2.1 The Spinor Representation 

The multiplet of states associated with topologically non-trivial solitons, on which the 
Clifford generators bi act, must contain a factor transforming according to the spinor 
representation of 0(2N). Since according to our discussion in subsection A. 1.1 (and 
according to [31]), the Dirac-Majorana equation (D.4) with a topologically non-trivial 
cr(x) has only a single zero energy bound state fo(x), there is only one such factor of 
the spinor representation associated with a given topologically non-trivial a(x). For 
example, CCGZ kinks, for which (D.4) has only a zero energy bound state, are pure 
isospinors[30]. 

It will be useful at this point to list a few mathematical facts concerning this 
0(2N) spinor representation factor: In any representation of (D.8) we may define the 
antihermitean 0(2N) generators as 

Mij = bibj - bjbi . (D.10) 

In addition, let us define the hermitean operator 1^ = 26162 • • • 62^, which commutes 
with all group generators M^, but anticommutes with all generators of the Clifford 
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algebra bi. Clearly, T\ = 1, and thus the possible eigenvalues of T 5 are ±1. There 
is only a single irreducible representation of the Clifford algebra (D.8), since T 5 does 
not commute with the bi, and any representation of the Clifford algebra (D.8) must 
contain both "left-handed" isospinors with T 5 = +1, and "right-handed" isospinors 
with r 5 = — 1. 

We would like to construct a unitary irreducible representation of this Clifford 
algebra on a positive norm Fock space. Thus, we have to group the Clifford generators 
bi into N pairs of Grassmannian creation and annihilation operators which act on that 
space. This we achieve e.g., by pairing the bi according to (D.3), namely 

B a = ^(b 2a -i + ib2a), a=l,...N, (D.ll) 

which satisfy the anticommutation relations 

{B a ,Bl} = 5 ab , (D.12) 

with all other anticommutators vanishing. An orthonormal basis for this space is then 
obtained by applying the creation operators B\ repeatedly on the Fock vacuum: 

|0>, flt|0>, BlBl\0), BlBt-.-BUo). (D.13) 

Thus, this irreducible representation of the Clifford algebra has dimension 2 N . 

However, since T 5 commutes with all the 0(2N) generators My, these 2 N states 
transform reducibly under the 0(2N) group; the left-handed and right-handed states 
transform independently. For example, if we choose the Fock vacuum such that 
r 5 |0) = +|0), then the right-handed states are those with an even number of B]fs 
acting on |0), and the left-handed states are those with an odd number of B^s acting 
on |0). For later reference, let us demonstrate this reducibility in a concrete basis of 
the algebra, in which the 0(2N) Cartan subalgebra is generated by the N generators 

M 12 , M 34 , ... M (2N _ 1){2N) . (D.14) 

From (D.ll) we see that M^ a -i){2a) — 2&2a-i&2a = — 2i[Bl,B a ]. Thus, the states 
(D.13) are simultaneous eigenstates of all the generators (D.14), and the correspond- 
ing eigenvalues are essentially the components of the weight vectors of the represen- 
tation. The step operators, on the other hand, when acting on any of these states, 
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change the number of Fock quanta of that state by or ±2 units (consider, e.g., 
M 13 = 2bib 3 = 2(Bi + B{)(B2 + B\)). In other words, the states with an even num- 
ber of B^s acting on |0) and the states with an odd number of B^s acting on |0) 
form two disjoint (and irreducible) invariant subspaces under the action of the 0(2N) 
generators. 

D.2.2 The Antisymmetric Tensor Representations 

Suppose now that (D.4), in a given solitonic background cr(x), has a pair of bound 
states at ±u n . What are the irreducible 0(2N) factors in the multiplet of states 
born by the soliton cr(x), which are acted upon by the operators a n ^ and aj^? To 
answer this question consider the following situation: Assume for simplicity that ±u n 
are the only bound states in (D.4). Then choose k distinct operators among the 2N 
£i's, and apply them to the vacuum state |0). The resulting state IljLi will 
contain bound as well as scattering states. Since we are interested only in static 
stable soliton states, we have to consider only those ai/s which correspond to bound 
states of (D.4). Thus, projecting out all scattering states we are left with the state 
(fn) k a n,ii ■ ■ ■ a n,i k 1^) (where {fn) k stands for the direct product of the k spinors). 
Clearly, we can form (2N)\/k\(2N — A;)! bound states in this way, all of which have 
common energy ku n , and thus form an energy multiplet. We readily identify this 
multiplet as the antisymmetric tensor of 0(2N) of rank k, since all a) n ^ anticommute. 
In this construction, we can obviously take 1 < k < 2N, and thus obtain all 0(2N) 
antisymmetric tensors. However, in the context of studying extremal static a(x) 
configurations in the GN model, the number v n of quanta trapped in the bound 
states ±u; n , is part of the definition of o~(x), in which case one must take k = v n . 

For pedagogical clarity, let us derive this conclusion in a more pedestrian way. 
To this end we have to invoke the notions of particles and antiparticles. In order 
to differentiate "particles" from "antiparticles" in our model we have to pair the 2N 
Majorana fields in (D.6) into N Dirac fields ip a . (The fermion number operator is 
of course / dx : SaV'lV'a : •) One possible pairing is according to (D.3), namely 
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ipa — ^2a-i + i&a, a = 1, . . . N 25 , which gives rise to the Dirac field operators 

i>a = 6a- 1 + l&a = fo(x)(b 2a -l + «&2a) + 

)t} . (D.15) 

w„>0 

Suppose that the Dirac-Majorana equation (D.4) has a pair of charge conjugate 
bound states at ±u n . Consider the (p n , h n ) many-body configuration of p n particles 
and h n holes defined following (3.20). The energy of this state is simply (p n + h n — 
N)ui n . It depends on the sum v n = p n + h n of numbers of particles and holes, and not 
on each one of them separately. Evidently, Fermi statistics implies that < v n < 2N. 
Thus, considering the set of all (p n , h n ) configurations with u n = p n + h n fixed, 
we obtain a multiplet of ££=0 C%C^ ~ k = = {2N)\/v n \(2N - v n )\ states with 
common energy {y n — N)uj n . This is, of course, precisely the result we obtained 
earlier directly from the 0(2N) symmetry, without distinguishing between particles 
and antiparticles. The many-body wave functions of states in this multiplet are 
Slater determinants of v n single particle orbitals, with 2N possible orbitals (i.e., with 
a flavor index running through 1, . . . 2N). Thus, we readily identify this multiplet as 
the antisymmetric tensor representation of rank v n of 0(2N). Since < u n < 2N, all 
possible 0(2N) antisymmetric tensor multiplets may occur. Thus, in particular, both 
the antisymmetric tensors of rank v n , and its dual tensor, of rank 2N — u n , appear 
in the spectrum. These two tensor multiplets have the same dimension and opposite 
energies ±(z/ n — N)u n . Note, however, that the extremal fermion bags, discussed 
in subsection 3.1, realize only half the antisymmetric tensors allowed by the 0(2N) 
symmetry, namely, only tensors of ranks in the range < u n < N (recall (3.26)). 

25 This is by no means the only possible pairing. There is clearly a one to one correspondence 
between permutations of the 2N fields £j and their pairings into N Dirac fermions, and thus 
(2iV)! ways to pair the 2N Majorana fields into N Dirac fields, namely, i\) a = £p(2 -i) + *£p(2a) = 

/o(&P(2a-l) + ib P(2a))+J2 Un >0 {fn( a n,P(2a-l) + *«n,P(2a)) + fn( a n,P(2a-l) - ia n ,P(2a)V} , where P 

is a permutation of 2N objects. Then clearly, the operator a„ p( 2a _i) + ia n .p(2a) ma y be thought 
of as annihilating a particle in the state f n (x,t), and the operator (a„ ; p(2o-i) ~ * a n,P(2a)) t ma Y 
be thought of as creating an antiparticle in the state f*(x,t). However, since the symmetry of the 
system in question is 0(2N) (rather than U(N)), all these pairings are physically equivalent. They 
simply correspond to different orientations of the corresponding U(N) subgroup in 0(2N). Thus, 
with no loss of generality, we choose to pair according to (D.15), and define the notions of particles 
and antiparticles (or holes) accordingly. 
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Referring, as usual, to the Dirac particles occupying the state at +oj n as fermions 
and to the Dirac holes in the state at — u n as antifermions, then the (p n , h n ) config- 
uration would carry fermion number = p n — h n . This is what we referred to as 
the valence fermion number N^ al in (3.43). Thus, states in our antisymmetric tensor 
multiplet carry different fermion numbers = p n — h n = v n — 2h n , which vary in 
the range -v n < ivf n) < v n (if u n < N), or -(2N-u n ) < ivj n) < 2N -v n (if v n > N). 
Thus, the spectrum in the antisymmetric multiplet of rank v n coincides with the 
7Vj n) spect rum in the dual antisymmetric multiplet of rank 2N — u n . 

Similarly, according to our explicit construction of the spinorial representation in 
the previous subsection, we see that the states in that representation carry different 
fermion numbers in the range < ]\fj spmoT ^ < JV; just count the number of B^s in a 
given state. (In subsection 3.2 we referred to this quantity as the valence fermion num- 
ber n .) However, this is an intolerable assignment of fermion numbers in a multiplet 
appearing in the spectrum of a quantum field theory with charge conjugation invari- 
ance, such as the GN model. One may say that the theory preserves its invariance 
under charge conjugation through the phenomenon of fermion number fractionaliza- 
tion [31, 33], which shifts, as we discussed in subsection 3.2 (recall (3.45)), the valence 
fermion numbers in the spinorial representation by —N/2, rendering the spectrum of 
fermion numbers in that representation symmetric: —N/2 < jyj spmor ' ) < N/2. More 
formally, we can identify the fermion number operator jvj spmor ' ) as the linear symmet- 
ric combination [37] 

• N N AT 

< Pin ° r) = J E M (2 _ 1)(2o) = £ B\B a - £ (D.16) 

* a=l a=l Z 

of the generators (D.14) of the Cartan subalgebra of 0(2N), where we used (D.12). 
The splitting of jvj spmor " ) into valence and fractional parts is manifest. 

D.3 Multiplet Dimensions from the Partition Function 

The dimensions of degenerate energy multiplets occurring in the quantization of 
(D.2) can be read off the partition function Z(/3) = Tr exp —f3H in a straightfor- 
ward manner. We consider here a topologically non-trivial cr(x), so that all possible 
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multiplets contribute. The Hamiltonian (D.9) is the sum of contributions of non- 
interacting Grassmannian oscillators, and the partition function is a product over the 
contributions of individual modes. The modes at frequency uj n contribute the factor 
Hi Tr exp —(3uj n a} n { a n ^ = (1 + exp —(3u n ) 2N to the partition function. The Fock space 
associated with the zero mode is 2 N dimensional, and thus contributes a factor 2^ to 
the partition function. 

Gathering all these facts together, we thus obtain the formal expression 

Z(P) = 2 N J] (l +e -^) 2N 

= 2 N n (e tw^t H ■ < D ' 17 > 

Thus, a state in which there are v n quanta at frequency uo n contributes a factor 
= (2N)\/u n \(2N— v n )\ to the overall degeneracy. This is, of course, the dimension 
of the antisymmetric tensor multiplet of rank v n . From the expansion of (D.17), 

Z{(3) = £ [2 N ]J C&) e~^i>^ 1 , (D.18) 

where the summation runs over sets of integers < Uk < 2N, we see that the 
eigenstates of the Hamiltonian (D.9) fall into all possible direct products of the 0(2N) 
antisymmetric tensors and the spinorial representation, which are built on the normal 
modes in the expansion (D.6) of £. 

It is somewhat amusing to obtain these results in the language of Dirac fields and 
the U(N) group. In this language it is more convenient to consider the Hamiltonian 
without normal ordering, i.e., to include the zero point energy (which we can subtract 
in the end). Thus, consider one of the iV-fold degenerate normal modes at frequency 
+u) n . If that state is empty the energy is — (i.e., the zero point energy), and if 
it is occupied the energy is Thus, these modes contribute a factor (exp ^ + 

exp— S^) N = (2 cosh to the partition function. The modes at — u n obviously 

contribute the same factor. Thus, the states at frequency ±u n contribute altogether 
a factor (2cosh^) 2iV = (1 + exp -f3uj n ) 2N expN[3u n . The TV-fold degenerate zero 
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mode contributes a factor (2 cosh/? • 0)^ = 2^, which is, of course, the dimension of 
the spinor representation. Multiplying all these factors together, and subtracting the 
overall divergent contribution of the zero point energies, we obtain (D.17). 
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Appendix E: Derivation of the Mass Formula for 

the CCGZ Kink 



In this Appendix we provide our own derivation of the mass formula M kink = 
(Eq. (3.33)) of the CCGZ kink. The kink mass, as all other soliton masses in the GN 
model, ought to be proportional to N. In addition, due to dimensional considerations, 
it has to be proportional to the dynamical mass m. Thus, on very general grounds, 

M kink = Nmc, (E.l) 

where c is a dimensionless pure number to be determined. In particular, c cannot 
depend on m. Therefore, 

1 dM kink ^ 



N dm 

M k i nk is given by the renormalized energy functional (2.19), evaluated at u k i nk {x) = 
m tanh mx. Note that due to the subtracted vacuum term in (2.19), that functional 
is also a function of uvac = m , and of course, also of the bare coupling g 2 . Thus, 

dMkink = 1 dx dakink(x) | dS da VAC | OS dg 2 
dm J 5a(x) dm dayAc dm dg 2 dm ' 

where all derivatives of £ are evaluated at cr(x) = cr kink (x). Now, a kink (x) and u VAC 
are extremal configurations. Thus, by definition 26 , (S£ /Sa(x)) akink ^ = dS/ddyAC — 
0. It follows from this and from (2.19) that 

dM kink mdg 2 1 d\ogg 2 



dm g 4 dm g 2 d log m ' 



;e.4) 



where we used / dx (<7 kink — m 2 ) = —2m. 

What is the meaning of taking the derivative of the bare coupling constant g 2 

with respect to the physically observable dynamical fermion mass ml The answer 

26 This simple argument can be verified by substituting the kink resolvents B and C (A. 41) into 
(3.10), and by using the fact that V c = a 2 + a' = m 2 = <Jy AC for the kink. Then, the sum of the first 
two terms on the right hand side of (E.3) produces a convergent integral in x, which is proportional 
to the bare gap equation (C.4), and thus vanishes. 
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is that in order to produce the RG invariant mass m at the large distance scale, we 

TV 

have to tune the bare coupling g 2 according to (2.5), namely, Ae N a' 2 w = m . The 
RG invariant mass m thus parametrizes the RG trajectory of the bare coupling g 2 as 
a function of the cutoff A. Therefore, (dg 2 /dm)5m measures the change in g 2 as we 
move to a neighboring trajectory, keeping A fixed. Thus, substituting <91ogg 2 /<91ogm 
from (2.5) into (E.4), followed by substituting the resulting dMkink/dm into (E.2), 
we finally obtain 

c=-, (E.5) 

71 

thereby proving (3.33). 
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Appendix F: An Alternative Proof of Eq. (4.12) 



In this Appendix we present an alternative proof of (4.12). The idea is to consider 
the behavior of ^ 

/(e 1 ,...,e„) = X>me n (f.i) 

n=l 

over the hyperplane 

t ra : Q 1 + ... + Q K = a + ^r, (F.2) 

where < r < K is an integer, and < a < |. This hyperplane corresponds, of 
course, to solitons with a fixed value of the maximal fermion number A^ a:E (.Dp arcnt ). 

We will now prove that . . . , ©x) attains its absolute minimum on S ra in the 

positive orthant, at the vertices of the intersection of S rQ and the hypercube [0, 
namely, the points 

= \ (<W + $nn 2 + ■■■ + S njlr ) + a5 nrir+1 , (n = 1, . . . , K) , (F.3) 

with all possible choices of r + 1 coordinates ii, . . . ,i r +i out of K. Once we have 
established that, the fact that the potentially stable solitons are given by (4.12) 
follows in a straightforward manner. 

By standard constrained extremum analysis one can show that the local extremum 
of /(©i, . . . , Qr) occurs at the symmetric point n = j^(a + |r),Vn, and that it 
is a maximum. Thus the minimum occurs at the boundary of the domain under 
consideration. By the same argument, any local extremum on the boundary is again 
a maximum. Thus, repeating this analysis, we conclude that the absolute minimum 
occurs at the vertices (F.3) of the intersection of S rQ and the hypercube [0, 
where 27 

min/i- =r + sino;. (F.4) 

An alternative elegant proof of this minimum behavior of (F.I), due to Raphael 

Yuster, goes as follows: Consider a sequence < ©i < . . . < ©x < |, subjected 

27p or r — o, (F.4) reduces to the well known inequality sin(^^ =1 0„) < X)^Li snl ®™ where 
Eli©n = a<|. 
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to (F.2). Assume that for some i, < 0, < i+ i < f . We will show that there 
exists another sequence of 0's, with the same sum, but with a lower sum of the sines. 
Thus, let 5 > be chosen such that 0; — 5 > and i+1 + S < f , i.e., < 5 < 
min{0j, ^ — 0j+i}. Modify the sequence under consideration by replacing 0« by ©« — 5 
and 0j+i by 0j+i + <5, keeping the other K — 2 terms unaltered. The new sequence thus 
obtained has the same sum as the original sequence, and thus defines another point 
on E m . We must show that D(5) = /(original sequence) — /(new sequence) > 0. 
Indeed, D(S) = sin©; - sin(0; - 5) + sin0 m - sin(0 m + 5). Clearly, D(0) = 0, 
and also D'(5) > in the relevant range of 5. Thus, D(5) increases monotonically 
with 5, and reaches its maximum at 5 max = min{0j, | — 0j + i}, where, depending on 
the initial condition at 5 — 0, either 0j — 5 max = or i+ i + 5 max = |. Thus, the 
sequence of 0's constrained to S rQ , which minimizes (F.l), cannot have more than 
one element in the interior of [0, |]. Thus, due to (F.2), the absolute minimum is the 
sequence in which the r largest 0's are |, one is a and the rest are zero, namely, 
the vertices (F.3). 

Thus, a parent soliton corresponding to a point in the interior of the intersection 
of E rQ , and the hypercube [0, , can decay into a final state with quantum numbers 
corresponding to the points (F.3), i.e., L = K and 9 n = 0^ in (4.3). In fact, such 
a parent soliton can also decay at least into the set of final states contained in small 
pockets above the vertices (F.3), which correspond to L = K and 8 n = + e in 
(4.3), with 

e < ^ Sn=i sin Q" ~ ( r + sin a ) 
K — r — 1 + cos a 

or into final states corresponding to L > K in (4.3), with 6>; = 0-^ for 1 < i < K 
and 9i = e for K + 1 < i < L, where 

Yln=i sin - (r + sin a) 
Sme< L^K • 

On the other hand, the parent soliton which corresponds to the vertices (F.3) has 
no open channel to decay through. Hence it is potentially stable. Indeed, if it could 
decay through a channel corresponding to 9 1 , . . . , 9l, then, from the requirement that 
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these parameters satisfy (4.3), we would have 

L 

^sin^ < r + sina. (F.5) 

i=\ 

Define the hyperplane 

S ra : e 1 + --- + 6 L = a + ^r. (F.6) 

From the analysis in this Appendix we know that the absolute minimum of J2f=i sin 9i 
over the intersection of £ ra and the hypercube [0, |] L is r + sin a. Thus, the points 
which satisfy the second inequality in (F.5) are bounded by Yh=\ @i < a + § r, in 
contradiction with the first inequality in (F.5). This completes our alternative proof 
of (4.12). 
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